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CHAPTER  I 
INTRODUCTION 


The  advancement  of  semiconductor  material  technology  has  made 
possible  the  development  of  new  semiconductor  devices  and  integrated 
circuits  (IC)  during  the  past  three  decades.  As  a result,  there  is  a 
remarkable  evolution  of  semiconductor  devices  from  the  early  discrete 
components  to  the  very  large  scale  integrated  circuits  (VLSI)  and  very 
high  speed  integrated  circuits  (VHSIC)  at  the  present  time. 

Shallow  impurity  dopants  are  intentionally  added  into  a 
semiconductor  via  thermal  diffusion  or  ion  implantation  process  in  order 
to  supply  conduction  electrons  or  holes  for  semiconductor  devices. 

Thus,  the  density  of  doped  impurities  will  significantly  affect  the 

electrical  characteristics  of  semiconductor  devices.  Occasionally,  a 

semiconductor  is  found  to  be  contaminated  with  unexpected  impurities 
which  cause  the  failure  of  device  operation.  These  unexpected 
impurities  are  brought  into  a semiconductor  device  during  the 
fabrication  processes  of  an  integrated  circuit,  where  the  processes 
become  so  complex  in  the  modern  semiconductor  technology  in  which 

numerous  processing  steps  were  used  during  each  IC  fabrication. 

In  the  design  and  fabrication  of  a high  speed,  low  power,  VLSI  or 
VHSIC,  it  is  necessary  to  increase  the  density  of  shallow  impurity 

dopants  in  order  to  reduce  the  size,  the  parasitic  capacitance,  and  the 
parasitic  resistance  of  active  devices  in  the  circuits.  In  addition. 
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even  a small  amount  of  contaminated  impurities  inside  a semiconductor 
will  significantly  affect  the  electrical  performance  of  active  devices 
and  reduce  the  IC  yield. 

The  main  objectives  of  this  research  are  (1)  to  develop  theoretical 
models  for  heavily  doped  semiconductor  devices  (e.g.,  pn  diodes, 
heterojunctions,  transistors,  FETS  . . . ) taking  into  account  the  heavy 
doping  effects  such  as  carrier  degeneracy,  band  gap  narrowing,  and 
dopant  density  dependent  dielectric  constant;  (2)  to  theoretically 
analyze  the  effects  of  heavy  doping  on  the  semiconductor  device 

characteristics,  and  (3)  to  perform  experimental  study  of  the  effects  of 
metallic  contamination  on  the  device  characteristics . 

A detailed  theoretical  study  of  the  effects  of  increasing  the 
density  of  shallow  impurity  dopants  on  semiconductor  devices  and  an 
experimental  study  of  the  metallic  contamination  on  device 

characteristics  have  been  carried  out  in  this  disseration;  these 
i ncl ude : 

(1)  Expressions  for  the  Debye  length,  built-in  potential  and 
depletion  layer  width  of  a heavily-doped  pn  junction  diode  as 
well  as  a heterojunction  diode  are  reformulated.  The  heavy 
doping  effects  such  as  carrier  degeneracy,  bandgap  narrowing, 
and  dopant  density-dependent  dielectric  constant  are 
accounted  for  by  using  the  expression  for  Fermi-energy  given 
by  [1],  and  bandgap  narrowing  model  proposed  by  [2],  as  well 
as  the  dopant  density  dependent  dielectric  constant  given  by 
[3]. 
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(2)  Theory  for  a heavily-doped  pn  junction  is  modified  by 
treating  it  as  an  anisotype  graded  heterojunction  and 
considering  the  bandgap  narrowing  effect.  This  effect  will 
reduce  the  bandgap  of  a heavily-doped  semiconductor  and 
result  in  the  bandgap  difference  between  p-  and  n-region. 
Thus,  the  theory  for  an  anisotype  graded  heterojunction  can 
be  modified  to  describe  the  electrical  behavior  of  a heavily- 
doped  pn  junction  assuming  a transition  region  from  one  side 
to  the  other.  This  transition  region  extends  beyond  the 
space-charge  region,  and  takes  place  over  a large  distance  as 
compared  to  interatomic  distances. 

(3)  Computer  analysis  of  the  effects  of  heavy  doping  on  various 
semiconductor  devices  such  as  pn  diode,  heterojunction , 
metal -semiconductor  field-effect  transistor  (MESFET),  and 
junction  field-effect  transistor  (JFET)  are  performed  using 
the  theoretical  models  presented  in  this  dissertation. 

(4)  The  effects  of  metallic  contamination  on  the  electrical 
characteristics  of  silicon  devices  are  experimentally 
studied.  The  experimental  tools  employed  in  this  study 
include  the  current  voltage  (I-V),  capacitance  voltage  (C-V), 
thermally  stimulated  capacitance  (TSCAP)  and  deep  level 
transient  spectroscoopy  (DLTS)  measurements.  From  these 
measurements,  one  can  determine  the  energy  levels,  density  of 
defects,  defect  profile,  and  the  capture  cross-section  for 
each  trap  level.  These  trap  levels  are  introduced  into  the 
forbidden  bandgap  of  semiconductor  by  the  contaminated 
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impurities.  Thus,  it  provides  us  with  means  of  identifying 
the  unknown  impurities.  In  addition,  we  can  study  the 
correlations  between  the  densities  of  metallic  contaminants 
and  the  electrical  characteristics  of  devices. 

Study  of  the  effects  of  heavy  doping  on  the  built-in  potential, 
Debye  length,  depletion  layer  width,  and  depletion  capacitance  in  a 
heavily-doped  abrupt  or  linearly  graded  pn  junction  diode  are  described 
in  Chapter  II.  The  effects  of  heavy  doping  on  the  Einstein  relationship 
for  both  majority  and  minority  carriers  in  a semiconductor  are 
theoretically  analyzed  in  Chapter  III.  The  theory  for  current-voltage 
relationship  of  a heavily-doped  pn  diode  is  proposed  in  Chapter  IV,  in 
which  the  diode  is  treated  as  a device  with  a linearly  varying  bandgap 
across  the  transition  region. 

In  Chapter  V,  the  effects  of  the  heavy  doping  on  the  built-in 
potential,  depletion  layer  width,  and  depletion  capacitance  in  a heavily 
doped  abrupt  or  graded  heterojunction  are  depicted.  Calculations  are 
also  performed  for  a AlGaAs-GaAs  heterojunction  diode  using  the 
expressions  derived  in  this  chapter. 

Chapter  VI  deals  with  the  effects  of  heavy  doping  on  current- 
voltage  characterisitcs  of  a heavily-doped  JFET  and  MESFET  device, 
computer  calculations  are  carried  out  for  a typical  GaAs  JFET  and  MESFET 

1 O O 

with  channel  doping  up  to  5 x 10iO  cm  . 

The  experimental  methods  employed  for  studying  the  effects  of 
metallic  contamination  on  the  device  performance  are  described  in 
Chapter  VII.  Chapter  VIII  presents  the  experimental  results.  A summary 
and  coclusions  are  given  in  Chapter  IX.  Suggestion  for  future  research 


are  also  included. 


CHAPTER  II 

EFFECTS  OF  HEAVY  DOPING  ON  THE  BUILT-IN  POTENTIAL, 

DEBYE  LENGTH,  DEPLETION  LAYER  WIDTH,  AND  DEPLETION 
CAPACITANCE  IN  A HEAVILY-DOPED  PN  JUNCTION  DIODE 

2.1  Introduction 

Since  heavy  doping  effects  play  an  important  role  in  the 
performance  of  semiconductor  devices  containing  heavily-doped  pn 
junctions,  many  experimental  and  theoretical  studies  [1-7]  have  been 
reported  on  determining  the  impurity  density-dependent  dielectric 
constant  [4],  bandgap  narrowing  [5-7]  and  analytical  modeling  of  these 
effects.  For  example,  Castner  et  al  . [4]  observed  a donor  density 
dependent  dielectric  constant  from  capacitance  measurement  on  heavily- 
doped  n-type  silicon.  Andrews  et  al  . [3]  proposed  an  analytical  model 
to  account  for  the  effect  of  impurity  density-dependent  dielectric 
constant  on  the  electric  field  and  charge  density  in  a graded  n+~p 
junction.  Dumke  [5]  determined  the  bandgap  narrowing  in  the  heavily- 
doped  n-type  silicon  from  luminescence  measurements,  and  confirmed  that 
the  density  of  states  in  the  shifted  bandedge  states  for  degenerate 
silicon  may  still  be  approximated  by  using  parabolic  density  of  states 
without  introducing  significant  errors.  However,  there  are  no  reports 
on  extension  of  the  non-degenerate  pn  junction  theory  to  account  for  the 
heavy  doping  effects  such  as  carrier  degeneracy,  impurity  density- 
dependent  dielectric  constant,  and  bandgap  narrowing  on  the  Debye 
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length,  built-in  potential,  the  depletion  layer  width  and  depletion 
capacitance  of  a heavily-doped  pn  junction  diode. 

The  purpose  of  this  chapter  is  to  analyze  the  effects  of  heavy 
doping  on  the  built-in  potential,  Debye  length,  depletion  layer  width, 
and  depletion  capacitance  in  a heavily-doped  abrupt  or  linearly  graded 
pn  junction  diode.  Calculations  were  performed  for  the  p++-n  abrupt 
junction  and  linearly  graded  junction  using  the  empirical  approximation 
for  the  reduced  Fermi -energy  given  by  [1],  the  empirical  model  for  the 
dopant  density-dependent  dielectric  constant  Ksn  given  by  [3],  and  the 
bandgap  narrowing  model  proposed  by  [2],  Details  of  all  these  models 
are  described  in  section  2.2.  Section  2.3  deals  with  theoretical 
analysis  and  reformulation  of  the  expressions  for  the  built-in 
potential,  Debye  length,  depletion  layer  width,  and  depletion 
capacitance  for  a heavily-doped  pn  junction  diode.  Results  and 
discussion  are  given  in  section  2.4.  Conclusions  are  given  in  section 
2.5. 

2.2  Models  for  the  heavy  doping  effects 

Analysis  of  the  heavy  doping  effects  on  diode  characteristics  of  a 
heavily-doped  pn  junction  was  made  possible  by  employing  the  results 
reported  in  [1-3]  for  the  heavy  doping  effects  such  as  carrier 
degeneracy,  dopant  density-dependent  dielectric  constant  and  bandgap 
narrowing.  Expressions  for  calculating  the  different  heavy  doping 

effects  are  described  as  follows: 


S 
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2.2.1  The  carrier  degeneracy 

The  effect  of  carrier  degeneracy  may  be  taken  into  account  by 
tractably  using  an  empirical  approximation  for  the  reduced  Fermi -energy , 
n,  given  by  [1].  Thus,  the  reduced  Fermi-energy  for  electrons  and  holes 
in  a heavily-doped  n-  and  p-type  semiconductor  may  be  approximated  by 


- 


F(Y„) 


(2.1) 


"p  - *n(Yp)  + F(Yp) 


where 


,1/2 


F<Vp>  ’ Yn,p  ' + 3'5354  ' Vn,p  ' + Ti*”'1'4 


(2.2) 

(2.3) 


and  Yn  = n/Nc,  Yp  = p/Nv,  n(p)  is  the  electron  (hole)  density;  Nc  (Nv) 
is  the  effective  density  of  conduction  (valence)  band  states.  Values  of 
F(Y)  vary  from  10-3  to  1 for  n-type  silicon,  and  2 x 10-3  to  1.8  for 
p-type  silicon  for  dopant  densities  in  the  range  of  1 x 10^  to  1 x 103(3 
cm"3.  Equation  (2.3)  may  generate  an  error  An  with  several  alterations 
of  sign  in  the  range  0 < n < 12,  but  with  peak  excursions  |Anmax|  < 10"4 
[8],  The  calculated  values  of  F ( Y)  for  electrons  and  holes  based  on 
Eq.  (2.3)  for  silicon  are  shown  in  Fig.  2.1.  Note  that  the  first  term 
on  the  r.h.s.  of  Eqs.  (2.1)  and  (2.2)  is  the  direct  result  of  the 
nondegenerate  statistics  which  relates  the  carrier  density  to  the 
reduced  Fermi  energy,  while  the  second  term,  F(Yn)  may  be  considered  as 
a correction  for  the  Fermi-Dirac  statistics  for  the  degenerate  case. 
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2.2.2  The  dopant  impurity  density-dependent  dielectric  constant 
Experimental  results  [4]  show  that  increase  in  the  static 
dielectric  constnat  for  Sb,  P,  and  As  donors  in  Si  was  associated  with 
increase  in  net  ionized  donor  density,  Nxn  = Nq  - N^.  The  results 
indicate  a possible  divergence  in  the  static  dielectric  constant,  Ksn  at 
a critical  concentration  NQ.  Typically,  NQ  ~ 4 x 10^8  cm"3  for 
phosphorus  and  ~ 7 x 1018  cm"3  for  arsenic.  For  example,  the  dielectric 
constant  for  N-type  Si  at  room  temperature  can  be  written  as  [3] 


K = 

11.7 

for  n < 

1016 

-3 

cm 

sn 

A exp  (BN)  , 

n r ' n v n ' 7 

for  n > 

1016 

-3 

cm 

By  using  Ksn  = Km  at  Nxn  = NQ,  and  Ksn  = 11.7  at  Nxn  = 1016  cm"3,  we 
obtain 

-K 

An  = 11.7  (km/11.7)  (2.5) 

Bn  = Kq  • 10“16  • *n(km/11.7)  (2.6) 

The  value  of  K , which  determines  the  rate  of  increase  in  K$n(n)  with 
dopant  density,  was  chosen  as  50.  Values  of  KQ  are  equal  to  0.002506 
and  0.001429  for  P and  As  doped  Si,  respectively.  The  model  given  by 
Eq.  (2.4)  is  sufficient  to  show  trends  and  its  consistency  with  the 
experimental  data  reported  by  [4],  in  which  values  of  K$n  are  varied 
from  11.7  to  72  for  phosphorus  doped  Si  and  11.7  to  33  for  As  doped 
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silicon  for  dopant  densities  in  the  range  of  1017  - 5 x 1018  cm-3, 
respectively. 

2.2.3  The  bandgap  narrowing 

Recently,  Berggren  and  Sernelius  [2]  have  proposed  a rigorous 
theoretical  model  for  bandgap  narrowing  calculations  in  heavily-doped1 
many  valley  semiconductor.  They  attributed  the  bandgap  narrowing  to  the 
free-carrier-impurity  interaction,  the  correlation  energy,  as  well  as 
the  random  distribution  of  impurity  ions. 

The  calculated  bandgap  narrowing  values  presented  in  [2]  were  found 
in  excellent  agreement  with  the  bandgap  narrowing  data  determined  from 
the  luminescence  measurements  by  Dumke  [5]  for  the  heavily-doped  Si. 
Thus,  the  bandgap  narrowing  values  reported  by  [2,  5]  will  be  used  to 
estimate  its  effects  on  the  diode  parameters  of  a heavily-doped  pn 
junction . 

We  represent  the  dopant  density  dependent  bangap  narrowing  for 
silicon  by  a five  degree  polynomial  formulation  for  dopant  density,  Np, 
greater  than  1037  cm-3  [2,  5]: 

AEg  = Aq  + Axx  + A2x2  + A3x3  + A4x4  + A5x5  (2.7) 

where  x = log(ND/1017) , AQ  = 17,  Aj-  = 15.483244,  A2  = 7.1771826,  A3  = 
0.06657874,  A4  = -1.1285211,  and  Ag  = 0.40151523. 

The  expressions  given  by  Eqs.  (2.1),  (2.2),  (2.4),  and  (2.7),  which 
account  for  the  heavy  doping  effects  described  above,  will  be  used  to 
compare  the  values  of  built-in  potential,  Debye  length,  depletion  layer 
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width,  and  depletion  capacitance  obtained  by  conventional 
(nondegenerate)  expressions  from  those  predicted  by  Eqs.  (2.1)  through 
(2.7)  for  the  heavily  doped  pn  junction. 

2.3  Derivation 

2.3.1  Built-in  potential 

The  built-in  potential,  4>bi  = <t>n  - <f>p  (<|>p  < 0),  is  the  total 
potential  change  from  the  quasi -neutral  p-type  ragion  to  the  quasi- 
neutral  n-type  region  (as  shown  in  Fig.  2.2);  it  is  a function  of  the 
dopant  density  in  each  region.  From  Fig.  2.2,  4>n  and  <j>p  are  the 
potential  at  the  quasi -neutral  edges  xn  and  -xp,  respectively,  and  may 
be  defined  by 


♦ = -(-  EN  - ♦*) 

Yn  vq  i n' 


<f>  E 


IP  * 
— E - + 4> 
q i P 


(2.8) 


(2.9) 


where  E^  (<  0)  and  E?  (<  0)  denote  the  intrinsic  Fermi  level  for  n-  and 
p-type  semiconductor,  respectively,  and 


★ 


- EN)  =*In 
c ' q n 


(2.10) 


★ 


(2.11) 


where  E^  and  E^  are  the  valence  and  conduction  band  edges  in  the  quasi- 
neutral p-  and  n-type  regions,  respectively.  Thus,  we  may  obtain  the 
modified  built-in  potential,  4>bi , for  a heavily-doped  pn  junction  as 
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and 


(2.12) 


(2.13) 


narrowing  in  the  n-  and  p-  type  material,  respectively. 

2.3.2  Debye  length 

In  order  to  investigate  the  validity  of  the  depletion  approximation 
for  a heavily-doped  abrupt  junction,  we  study  the  effects  of  heavy 
doping  on  the  extrinsic  Debye  length  by  using  the  standard  procedures 
described  elsewhere  [9]. 

If  we  consider  a uniformly  doped  p++-n+  junction  with  the  space 
charge  region  edges  at  -Xp  and  xn,  assuming  only  small  variations  of 


potential  from  <j>  near  x = xn  and  neglecting  the  minority  carrier 


* 


concentration,  p,  with  <f>'  = <j> 


n 


<f>,  we  may  then  write  the  Poisson's 


equation  in  the  form 


l/2vkTVHn 


(2.14) 
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where  is  the  Fermi-Dirac  integral  of  order  one-half,  and  esn  (=  Ksn 

* eQ)  is  the  permitivity  of  the  seimconductor;  eQ  is  the  permitivity  in 
vacuum,  and  Ksn  is  the  dopant  density-dependent  dielectric  constant 

given  by  Eq.  (2.4).  We  may  expand  ^ \j^  C^y  ( 4>n  - V)]  in  Eq.  (2.14)  in 
a Taylor  series,  and  retain  only  the  first  two  terms,  which  yields 


d2<j> 

dx2 


_a_ 

e 

sn 


F-l/2(kT  V 
kT  Fl/2(kT  *n} 


"DH 


(2.15) 


where  LDH  is  the  modified  extrinsic  Debye  length  for  a heavily-doped 
abrupt  pn  junction,  and  is  given  by 


, r£snkT  F1/2(V  n 
LDH  = F , , Jj 


1/2 


D 


- 1/ 2 v n' 


LD^F 


1/2^ 

-1/2^  V 


1/2 


■] 


(2.16) 


where  LD  is  the  Debye  length  for  the  nondegenerate  case.  It  is  noted 
that  Eq.  (2.16)  will  reduce  to  the  coventional  expression  as  dopant 
density  approaching  nondegenerate  regime.  Landsberg  [13]  has  also 
reported  similar  results,  but  without  including  the  effect  of  dopant 
impurity  density-dependent  dielectric  constant.  Thus,  it  is  obvious 
that  both  dopant  density-dependent  dielectric  constant  and  carrier 
degeneracy  will  have  significant  effects  on  the  extrinsic  Debye  length. 
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2.3.3  Depletion  layer  width  and  depletion  capacitance 
The  depletion  layer  width  and  depletion  capacitance  are  two 
important  parameters  for  the  design  and  modeling  of  a semiconductor 
device.  Since  the  dielectric  constant  is  a function  of  dopant  impurity 
density,  it  will  obviously  affect  the  depletion  layer  width  and 
depletion  capacitance  of  a heavily-doped  abrupt  junction,  as  well  as  a 
linearly  graded  junction,  and  result  in  a change  of  depletion  layer 
width  and  depletion  capacitance  from  that  given  by  the  nondegenerate 
case.  We  shall  first  deal  with  a uniformly  doped  p++-n+  abrupt  junction 
(as  shown  in  Fig.  2.3).  If  e$n  (=  K$n  • eQ)  and  e$p  (=  Ksp  * eQ)  are 
the  permitivity  in  the  n-  and  p-type  regions,  respectively,  and  Ksn  and 
K are  a function  of  dopant  impurity  density,  then  we  may  write 

b p 


E(x)  = (x[J  - x)  , 0 < x < x 


H 


n 


sn 


-xH  < x < 0 
P 


(2.17) 


sp 


Thus,  the  potential  is  given  by 


o sn 


(2.18) 


H H 

for  the  abrupt  junction,  with  • xp  = • xn,  the  depletion  layer 

width  for  a heavily-doped  abrupt  junction  can  be  expressed  by 
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x , = 


H , H 
x + x„ 
n p 


2e 

= [•  ° 


*bi 


q N. 


1/2  N.  1/2  Nn  1/2 

i • t<r>  + <nt>  ] 

inD  A 


(2.19) 


sn 


sp 


Since  the  depletion  layer  width  is  dominated  by  the  lightly-doped 
region,  we  may  write  the  depletion  layer  width  for  the  heavily-doped 
abrupt  junction  as  a function  of  applied  voltage  by  including  the  effect 

of  two  majority-carrier  distribution  tails  [10]  and  replacing  the 

u u 2k  T 

modified  built-in  potential  <f>b^  in  Eq.  (2.19)  by  ( 4>bl-  - Vg — ) so 

that 


= [- 


2e 


u 


bi 


V - 
a 


2k  T\  1/2 


(NnTK 

v A sn 


Nn/k  ) 
D'  sp' 


■] 


NA 

[(/) 

ind 


1/2 


ND 

(/) 

ina 


1/2 


(2.20) 


The  depletion  capacitance  per  unit  area  for  a heavily-doped  abrupt 
junction  is  given  by 


C .4 2_  « ni/zni/2  . (. 

jA  dV  D A 

3 < 


1/2 


<VKsn  + VKsp>  * Ki 


V - — ) 
a q 


) 


£ / X 

snp'  xd 


(2.21) 


where 


£.(Nfl  + 


'snP  " (VKsn  + 


(2.22) 
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which  shows  that  dielectric  constant  in  the  heavily-doped  n-  and  p-type 

regions  depends  on  dopant  impurity  density  instead  of  a constant  es  = 

11.7.  The  depletion  capacitance  for  the  heavily  doped  case  may  be 

obtained  from  Eq.  (2.21),  which  is  also  dopant  density  dependent. 

Similar  calculations  may  be  performed  for  the  case  of  a heavily- 

doped  linearly  graded  junction  with  Nq  - = ax,  where  "a"  is  the  net 

doping  concentration  gradient.  Detailed  mathematical  derivation  of  the 

H H 

depletion  layer  width  and  the  built-in  potential,  <|»^. , in  a heavily- 
doped  1 i nearly-graded  junction  taking  into  account  the  heavy  doping 
effects  are  described  in  the  Appendix.  It  can  be  shown  that  from  Eq. 
(A— 10 ) 


*bi 


<j>  - 4> 

yn  Yp 


- 


J 1 r 1 / H .2 

1 a ar  ) + 


L.  H 

H >2  , 1 rnxdL 


2e  1 A L 4C  VAdL' 
o n 


dL 


? Ir  H 

1 r 1 , H / , 1 "rpxdL  H 

'A  4C  r2  e * xdL 

P P cp 


lp  H 

—j  • (e^"XdL  - 1)] 

cJ 

n 

Jr  H 

% (e^L  - 1)]} 

S 


(2.23) 


where  a is  the  gradient  of  net  doping  concentration,  An  (Ap)  is  the 
constant  of  K$n  (K$p),  repsectively , and  Cn  = -Bna,  Cp  = Bpa  are  the 
parameters  characterized  only  by  junction  properties  such  as  the  dopant 
impurity  density  and  the  net  dopant  concentration  gradient,  a.  It  is 
noted  that  the  depletion  layer  width,  x^J.  , may  be  expressed  as  a 
’function  of  applied  voltage  with  eg  = 11.7  (i.e.,  without  taking  into 
account  the  effect  of  dopant  density-dependent  dielectric  constant)  by 


(2.24) 
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PL 


+ x 


nL 
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/ , H 

t ( 't’u  • 

s v 


qa 


1/3 

•] 


Thus,  the  depletion  capacitance  for  a heavily-doped  linearly  graded  pn 
junction  with  e$  = 11.7  is  given  by  [11] 


= [- 


qesa 


12(  <p  - 

g 


1/3 

■] 


(2.25) 


where  the  gradient  voltage  [12]  for  the  heavily-doped  linearly  graded 
junction  is  given  by 


2 A Tx  r a £skT/q  _ 

3 lq  ; „P  x 3/ 2 


] 


(2.26) 


8q(n.  n.  )' 
le  ie' 


and  the  effective  intrinsic  density  nie  and  n^e  in  the  n-  and  p-regions 
are  given  respectively  by: 


n^e  = [n?  exp  ( AEg/kT)F  1/2 ^ nn ^ exp  ("V^2  (2.27) 

n^e  = [n?  exp  (AEg/kT)F1/2(np)  exP  (_Tip)]1/2  (2.28) 


where  ni  is  the  intrinsic  carrier  denisty. 

In  order  to  determine  the  depletion  layer  width  in  a heavily-doped 

linearly  graded  junction,  we  need  to  choose  a suitable  trial  value 

of  x"  . and  then  use  iteration  procedures  on  Eqs.  (2.12),  (2.23)  or 
dL 

(2.24)  with  a known  net  dopant  concentration  gradient,  a.  The  exact 

U 

value  of  x.,  can  then  be  obtained. 
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2.4  Results  and  Discussion 

Calculations  of  the  effects  of  heavy  doping  on  the  built-in 
potential,  Debye  length,  depletion  layer  width,  and  depletion 
capacitance  for  a heavily-doped  pn  junction  were  carried  out  by  using 
the  expressions  derived  in  section  2.1.  Figure  2.4  shows  the  calculated 
built-in  potential  vs.  background  doping  density,  Ng,  for  a p++-n  (or 
n++-p)  abrupt  junction.  The  doping  density  of  p++  (or  n++)  side  is 
assumed  to  be  1 x 102^  cnT^,  and  Ng  is  the  impurity  concentration  on  the 
lightly-doped  side.  From  Fig.  2.4,  it  is  shown  that  the  modified  built- 
in  potential,  xj\ , calculated  from  Eq.  (2.12),  is  about  50  to  75  mV 
smaller  than  those  calculated  from  the  conventional  expression.  Thus, 
the  heavy  doping  effects  such  as  carrier  degeneracy  and  bandgap 
narrowing  have  a significant  influence  on  the  value  of  built-in 
potential  of  a heavily  doped  abrupt  pn  junction. 

The  result  for  the  Debye  length,  depletion  layer  width  and 
depletion  capacitance  of  a p++-n  abrupt  junction  are  shown  in  Fig.  2.5, 
2.6  and  2.7,  respectively,  in  which  the  doping  density  of  the  p++  side 

on  o 

is  assumed  to  be  1 x 10cu  cm  . From  Fig.  2.5,  it  is  clearly  shown  that 
the  effects  of  carrier  degeneracy  and  dopant  density-dependent 
dielectric  constant  have  considerable  influence  on  the  Debye  length. 
The  effect  of  dopant  density-dependent  dielectric  constant  on  the  Debye 
length  is  clearly  shown  in  Table  2.1.  Figure  2.6  and  2.7  show  the 
depletion  layer  width  and  depletion  capacitance  calculated  from  Eqs . 
(2.20)  and  (2.21),  with  and  without  including  the  effect  of  dopant 
density-dependent  dielectric  constant.  The  results  shows  that  the 
effects  of  carrier  degeneracy  and  bandgap  narrowing  have  little 
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influence  on  the  depletion  layer  width  and  depletion  capacitance.  On 
the  other  hand,  the  result  shows  that  the  dopant  density  dependent 
dielectric  constant  is  the  only  heavy  doping  effect  which  will  have  a 
strong  infuence  on  these  two  parameters.  Table  2.1  shows  the  calculated 
results  of  LD  and  xn  at  Ng  = 5 x 1018  cm'3.  Thus,  the  depletion 
approximation  is  still  valid  up  to  this  dopant  density. 

Fig.  2.8  shows  the  gradient  voltage  calculated  from  Eq.  (2.26)  for 
a heavily-doped  1 inearly-graded  junction,  where  dielectric  constant 
e$  = 11.7  is  used.  Thus,  the  change  in  the  gradient  voltage  shown  in 
Fig.  2.8  may  be  attributed  to  the  effects  of  bandgap  narrowing  and 
carrier  degeneracy.  The  results  of  depletion  layer  width  and  depletion 
capacitance  for  a heavily-doped  linearly-graded  junction  calculated  from 
Eqs . (2.24)  and  (2.25)  are  shown  in  Fig.  2.9,  where  dielectric  constant 
e$  = 11.7  is  used.  From  Fig.  2.9  it  is  clearly  shown  that  carrier 
degeneracy  and  bandgap  narrowing  have  little  effect  on  the  depletion 
layer  width  and  depletion  capacitance.  The  results  show  little 
difference  from  those  predicted  by  the  nondegenerate  expressions. 

From  the  results  shown  above,  it  is  clear  that  the  doping  density- 
dependent  dielectric  constant  is  important  in  the  heavy  doping  regime. 
For  an  isotropic  semiconductor,  the  dielectric  constant  depends  on  the 
properties  of  atoms  and  the  manner  in  which  the  atoms  are  assembled  to 
form  a crystal  (through  the  effect  of  local  E-field).  Clausisus- 
Mossotti  relation  [4]  can  relate  the  dielectric  constant  to  electronic 
polarizability,  which  is  resulted  from  the  displacement  of  electron 
shell  (or  cloud)  relative  to  the  atom  nucleus.  Since  the  high  extrinsic 
impurity  concentration  and  the  nonuniform  distribution  of  these 
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impurities  will  probably  affect  the  local  E-field  and  alter  the 
displacement  between  nucleus  and  electron  shell,  the  heavily-doped 
semiconductor  dielectric  constant  is  also  expected  to  be  changed. 
Further  experimental  investigation  on  the  dielectric  constant  anomaly  is 
needed . 

2.5  Summary  and  Conclusions 

A theoretical  analysis  of  the  heavy  doping  effects  on  the  built-in 
potential,  depletion  layer  width,  Debye  length,  and  depletion 

capacitance  for  a heavily-doped  pn  junction  has  been  presented  in  this 
chapter.  Numerical  calculations  of  the  heavy  doping  effects  such  as 
carrier  degeneracy,  dopant  density-dependent  dielectric  constant,  and 
bandgap  narrowing  were  carried  out  by  employing  the  empirical  models 
described  in  section  2.2.  The  results  of  present  studies  show  that 
1)  bandgap  narrowing  and  carrier  degeneracy  have  important  effects  on 
the  junction  built-in  potential;  2)  carrier  degeneracy  and  dopant 
impurity  density-dependent  dielectric  constant  have  a significant 
influence  on  the  Debye  length  for  the  abrupt  junction  case,  and  3)  the 
dopant  impurity  density-dependent  dielectric  constant  is  the  only  heavy 
doping  effect  which  will  have  a strong  influence  on  the  depletion  layer 
width  and  depletion  capacitance.  The  present  results  showed  that  the 
new  calculated  values  for  the  Debye  length,  built-in  potential,  and 
depletion  capacitance  may  differ  significantly  from  those  predicted  by 
the  nondegenerate  theories. 
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TA8LE  2.1 


A comparison  of  the  calculated  diode  parameters  for  a heavily-doped 
p++-n  diode  with  and  without  including  the  heavy  doping  effects. 


*Physical  Nondegnerate  Degenerate  Degenerate  Degenerate 
Parameters  e,.  = constant  e = constant  As-doped  Si  P-doped  Si 

es(ND)  es(ND) 


LD(A) 

18 

25 

45 

63 

Xd(A) 

164 

160 

280 

420 

xn(A) 

156 

152 

267 

400 

★ 


p++  = 1.0  x 1020  cm-2,  n - 5 x 10^®  cm"2. 


FUNCTION  F(Y) 
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Fiqure  2.1  Values  of  function  F(Y)  vs  hole  or  electron  density 
in  the  heavily-doped  Si  at  T = 300  k.  Where  F(Y) 
is  calculated  from  Eq.(2.3).  Solid  line  is  values 
for  electrons,  dash  line  for  holes. 
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Figure  2.2  Energy  band  diagram  for  a heavily-doped  p -n  junction 
diode. 
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Figure  2.3 


Schematic  diagram  for  a heavily-doped  p -n  junction 
diode  used  in  this  study. 


Built-in  potential  (volt) 
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Figure  2.4  Built-in  potential  vs.  background  doping  density,  N„,  for 
a heavily-doped  silicon  abrupt  pn  junction  at  T = 300  K. 

Whefe,  the  doping  dnesity  of  p++(  or  N++  ) side  is  assumed 
to  be  1 x 10^u  cm-3,  and  N„  is  for  the  dopant  impurity 
concentration  in  the  lightTy-doped  side.  Solid  line  is 
the  predicted  values  by  the  conventional  theory  for  a p++- 
n junction,  dash  line  and  dash-dot  line  are  calculated  by 
using  Eq.  (2.12)  for  the  n++-p  and  p++-n  abrupt  junction 
d iode. 


Debye  Length  (A) 
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Figure  2.5  The  Debye  length  vs.  background  dopant  density,  N„,  for 
a heavily-doped  p++-n+  Si  abrupt  junction  at  T = 300  K. 
Solid  line  is  the  conventionally  predicted  values,  dash 
line  is  calculated  from  Eq.  (2.16),  using  dielectric 
constant  ( = 11.7  ),  while  dot-dash  line  and  double  dot- 
dash  line  are  calculated  from  Eq.  (2.16)  using  dopant 
impurity  density-dependent  dielectric  constant  for  the 
As-  and  P-doped  Si . 


Depletion  Layer  Width  (j*m) 
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Background  Doping  Density  (N  B ) (cm-3) 


Figure  2.6  Depletion  layer  width+^s .+background  dopant  density,  N , 
for  a heavily-doped  p -n  Si  abrupt  junction  at  Va  = 8 
and  T = 300  K.  Solid  line  is  the  conventional  predicted 
values;  dash  line  is  calculated  from  Eq.  (2.20)  with  di- 
electric constant  ( = 11.7  ).  While,  dash-dot  and  double 
dot-dash  lines  are  obtained  from  Eq.  (2.20)  using  dopant 
impurity-density-dependent  dielectric  constant  for  the 
As-  and  P-doped  Si,  respectively. 
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Figure  2.7  Depletion  capacitance  vs.  background  dopant  density,  NB, 
for  a heavily-doped  p++-n+  Si  abrupt  junction  at  Va  = 0 
and  T = 300  K.  Solid  line  is  the  conventionally  predicted 
values;  dash  line  is  calculated  from  Eq.  (2.21)  with  di- 
electric constant  equql  to  11.7.  While  dash-dot  and 
double  dot-dash  lines  are  obtained  from  Eq . (2.21)  using 
dopant  impurity  density-dependent  dielectric  constant 
for  the  As-  and  P-doped  Si,  respectively. 
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Figure  2.8  Graident  voltage  vs.  net  impurity  concentration  gradient, 
a,  for  a heavily-doped  Si  linearly  graded  junction  at  Va= 
0 and  T = 300  K.  Solid  line  is  the  conventionally  predi- 
cted values,  dash  line  is  calculated  from  Eq.  (2.26)  with 
dielectric  constant  equql  to  11.7. 
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Figure  2.9  Depletion  layer  width  and  depletion  capacitance  vs. 

net  impurity  concentration  graident,  a,  for  a heavily 
-doped  linearly  graded  Si  pn  junction  diode  at  Vg  = 0 
and  T = 300  K.  Solid  line  and  dash  lines  are  the  conven- 
tionally predicted  depletion  layer  width  and  depletion 
capacitance,  respectively;  dash-dot  line  and  double  dot- 
dash  lines  are  the  values  of  depletion  layer  width  and 
depletion  capacitance  calculated  from  Eqs.  (2.24)  and 
(2.25)  using  dielectric  constant  equal  to  11.7. 
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EFFECTS  OF  HEAVY  DOPING  ON  THE  EINSTEIN  RELATION 


3.1  Introduction 

Studies  of  the  effect  of  carrier  degeneracy  on  Einstein  relation 
for  the  majority  carriers  [13-15],  as  well  as  the  effect  of  bandgap 
narrowing  on  the  effective  intrinsic  density  [16]  and  the  minority 
carrier  diffusivity  [17]  in  the  heavily-doped  silicon  have  been  reported 
by  many  investigators.  The  reason  for  this  may  be  attributed  to  the 
poor  agreement  between  the  experimental  data  and  the  theoretical 
prediction  on  the  transport  parameters  in  the  heavily  doped  regime,  and 
lack  of  understanding  of  the  physics  underlying  the  basic  mechanisms  of 
heavy  doping  effects  in  the  degenerate  semiconductors.  Most  of  the 
previous  research  efforts,  however,  have  been  focused  on  heavy  doping 
effects  in  the  degenerate  silicon  due  to  its  importance  in  the  IC 
appl ications . 

Measurements  of  the  minority  hole  diffusivity  [18-22]  and  the  drift 
mobility  [23]  have  been  reported  recently  for  heavily-doped  n-type 
silicon.  Although  large  discrepancies  still  exist  among  these  reported 
values  for  the  minority  hole  diffusivity,  the  results  indicate  that 
values  of  the  minority  hole  diffusivity  tend  to  saturate  for  dopant 
density  exceeding  1019  cm-3  as  is  shown  in  Fig.  3.1. 

Determination  of  the  minority  hole  diffusivity  in  the  heavily-doped 
n-type  silicon  is  usually  carried  out  by  using  the  Einstein  relation 
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(i.e.,  Dp/yp  = kT/q)  assuming  that  the  majority  and  minority  hole 
mobilities  are  equal  in  the  heavy  doping  regime  [16,24].  However, 
recent  theoretical  calculations  [25],  and  experimental  observations 
[22,23]  show  that  the  values  of  minority  hole  mobility  are  considerably 
lower  than  that  of  majority  hole  mobility  in  the  p-type  silicon  for  the 
corresponding  dopant  density.  In  addition,  it  is  very  important  to 
investigate  how  the  effect  of  bandgap  narrowing  will  affect  the  Einstein 
relation  for  majority  carriers  and  minority  carriers. 

3.2  Theoretical  analysis 

3.2.1  Einstein  relation  for  majority  carriers 

The  effects  of  carrier  degeneracy  on  Einstein  relation  for  majority 
carriers  in  a degenerate  semiconductor  have  been  studied  by  employing 
parabolic  density  of  states  and  rigid  band  approximation.  These  studies 
proved  that  the  generalized  Einstein  relation  for  an  n-type  degenerate 
semiconductor  may  be  written  as 

-p  = (^-)[n/(dn/dn  )] 
yn  q 

• <r)[Fi/2(r,c)/F-i/2(nc)]  (3-!) 

Fig.  3.2  and  3.3  show  the  Dp/yp  and  Dp/yn  ratio  vs.  carrier  density  for 
a p-  and  n-type  degenerate  silicon,  respectively,  calculated  from  Eq. 
(3.1)  for  temperatures  from  100-400  K.  However,  the  electron  density  in 
a degenerate  semiconductor  can  be  expressed  by 
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n NcF1/2*"c)  Nc  r( 3/2) 


ill 

1 + exp(e  - 


de 


= N 


1 

c r ( 3/  2 ) 


/ 


(e‘  - ec) 


1/2 


de 


ec  1 + exp[(e 1 - ec)  - nc] 


(3.2) 


where  e = E/kT  is  the  reduced  energy,  ec  = Ec/k T is  the  reduced 

conduction  band  edge,  and  for  a heavily-doped  semiconductor,  Ec  = - 

AE  (N).  AEa(N)  is  the  bandgap  narrowing  and  depends  on  the  carrier 

y y 

density.  Since  ec  is  an  implicit  function  of  reduced  Fermi  energy,  nc, 
the  mathematical  derivation 


cF-1/2^c 


(3.3) 


needs  to  be  carefully  examined  when  the  bandgap  narrowing  effect  is 
considered. 

By  using  Leibnitz's  rule  for  differentiation  of  integrals,  we  may 

write 


dn 

dn. 


/ 


ec(n)  dnc  [1  + exp( (e 1 - ec)  - nc)] 

1 


d { [e‘  “ ec(n)] 


1/2 


-}  • de  1 


N . 


ec(n)  c 1 + exp(e ' - ec(n)  - nc) 

r , , N -,-1/2  dec(n)r  1 

- [e1  - ec(n)]  ^— [ 


c 1 + exp(e‘  - e - 


(e*  - ec(n)) 


-]}  * de' 


1/2 


nc] 


(3.4) 
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After  employing  integration  by  parts  and  steps  of  manipulation,  the 
integration  of  first  term  of  Eq.  (3.4)  can  be  expressed  by 


/ 

ec(n) 


d\  1 + exp(e 

oo 

- / Cl  + 

ec(n) 


— ] • (e-  - er(n))1/2  • de' 

- ec(n)  ' \] 

de  (n)  i(e'  - e (n))"1/2 

c 3 . C - 9 ] . de1 

anc  1 + exp(e 1 - ec(n)  - nc) 


(3.5) 


Substituting  Eq.  (3.5)  into  (3.4),  we  obtain 


dn 


c 1 


dn.  “ 7(3/2) 


2 -L  1 


-1/2 


■c  - 0 * + exp(e  - \] 


de 


= NcF-l/2<\> 


(3.6) 


This  proves  that  the  generalized  Einstein  relation  for  majority 
carriers  in  a heavily-doped  semiconductor  is  still  valid  if  the 
parabolic  density  of  states  and  rigid  band  approximation  were  assumed. 


3.2.2  Einstein  relation  for  minority  carriers 

It  is  instructive  to  study  the  effect  of  bandgap  narrowing  on 
Einstein  relation  for  minority  carriers.  The  PN  products  for  an  n-type 
degenerate  semiconductor  is  given  by 
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pn 


ru  • exp( AEg/kT) 


F1/2(nc)  • exp(-nc) 


(3.7) 


where  is  the  intrinsic  carrier  density;  nc  = (Ep  - Ec )/kT  is  the 
reduced  Fermi  energy;  n^e  is  the  effective  intrinsic  carrier  density; 
and  the  bandgap  narrowing  AEg  is  defined  by 

AE  = E1  - E (n)  (3.8) 

g g g 

where  Eg(n)  is  the  perturbed  bandgap.  With  this  definition,  AEg  is  a 
positive  quantity,  which  varies  with  n.  Taking  the  partial  derivative 
on  both  sides  of  Eq.  (3.7)  with  respect  to  nc  and  dividing  by  pn  yields 

(^-)/P  = (dnie/dV/nie  " (dn/dnc)/n  (3-9) 

c 

By  substituting  Eq.  (3.6)  into  Eq.  (3.9),  we  obtain 


$-)/p  = -C1  - (d4E  /dnc)(pr) 
c 3 


(3.10) 


The  generalized  Einstein  relation  for  minority  holes  is  given  by 


u 


P 


(— )[p/(dp/dnv)] 


(3.11) 


Using  the  chain  rule,  we  can  write 
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(3.12) 


Finally,  from  Eq.  (3.10)  to  Eq.  (3.12),  we  obtain 


= kT 


(3.13) 


Thus,  the  Einstein  relation  for  minority  carriers  in  a heavily-doped 
semiconductor  is  still  valid,  if  the  parabolic  density  of  states  and 
rigid  band  approximation  are  assumed. 

3.3  Discussion 

From  the  above  theoretical  analysis,  it  is  clearly  shown  that  the 
generalized  Einstein  relation  for  majroity  carriers  and  minority 
carriers  in  a heavily-doped  semiconductor  are  still  valid  when  the 
parabolic  density  of  states  and  rigid  band  approximations  are  assumed. 
Physically,  with  any  gradient  of  majority  carriers  concentration  there 
will  also  be  a gradient  of  the  bandgap  which  gives  rise  to  an  effective 
E-field.  This  E-field  must  be  taken  into  account  in  representing  the 
drift  current  of  the  minority  holes.  When  this  effective  field  is 
included,  the  Einstein  relation  for  the  minority  carriers  is  still  given 
by  the  usual  expression,  Dp  = (kT/q)Up.  The  reasons  for  the  difference 
between  the  diffusion  constants,  say,  holes  in  a heavily-doped  p-type 
semiconductor  and  as  a minority  carrier  in  a heavily-doped  n-type 
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semiconductor  may  be  due  to  a combination  of  1)  the  difference  in  the 
distribution  that  holes  would  have  in  both  types  of  semiconductor  and 
2)  the  addition  of  electron-hole  scattering  for  holes  as  minority 
carri ers . 

It  is  noted  that  Eq.  (3.1)  is  not  valid  when  considering  the 
impurity  band  broading  and  band-tail  states  [24,  25]  instead  of 

employing  the  parabolid  density  of  states  and  rigid  band 
approximation.  Additional  terms  should  be  added  in  Eq.  (3.1),  and 
problems  will  become  more  complicated.  Further  research  in  this  area  is 
needed. 
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Figure  3.1  Minority  hole  diffusivity  vs.  dopant  density  in  n-type  Si. 

Solid  dots  are  the  data  given  by  [18],  open  circles  by 
[19],  solid  squares  by  [20],  and  open  squares  by  [21]. 


/ Ap  tVo,t) 


38 


Figure  3.2  Majority  hole  diffusi vity/mobil ity  ratio  vs.  hole 

density  with  temperature  as  a parameter  for  p-type  Si. 
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Electron  Density  ( cm-3) 


Figure  3.3  Majority  electron  diffusivity/mobil ity  ratio  vs.  elec- 
tron density  with  temperature  as  a parameter  for  n-type 
sil icon. 


CHAPTER  IV 

THEORY  OF  ELECTRICAL  CHARACTERISTICS  OF  A HEAVILY-DOPED  PN  DIODE 


4.1  Introduction 

Bandgap  narrowing  is  one  of  the  fundamental  heavy-doping  mechanisms 
that  yields  a large  mi nority-carrier  density,  and  consequently  a large 
minority-carrier  emitter  current.  Electrical  measurements  in  heavily- 
doped  silicon  material  and  devices  indicate  that  the  forbidden  energy 
gap  may  be  smaller  than  the  intrinsic  energy  gap  of  pure  silicon.  These 
experimental  measurements  are  categorized  by  [26]  as:  1)  measurements 

of  the  temperature  dependence  of  the  current  gain  in  a silicon  bipolar 
transistor  [27,28];  2)  measurements  of  the  collector  current  and  its 
temperature  dependence  in  silicon  bipolar  transistors  [16,29]; 
3)  measurements  of  the  temperature  dependence  of  the  photoresponse  of 
silicon  pn  junctions  and  transistors  [30];  4)  measurements  of  the 

temperature  dependence  of  the  emitter  current  in  silicon  solar  cells  (in 
the  dark  condition)  and  bipolar  transistors  [31];  5)  measurements  of  the 
emitter  current  in  different  pn  jucntion  structures  [20];  and 
6)  measurements  of  the  collector  current  in  different  pnp  bipolar 
devices  with  different  base  width  [32]. 

Rather  than  providing  a direct  measurement  of  the  bandgap,  the 

2 

above  measurements  yield  an  effective  intrinsic  carrier  density,  n.j  , by 
employing  the  conventional  current-voltage  expression.  However,  a 
question  arises  as  to  whether  the  heavy  doping  effects  such  as  bandgap 
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narrowing,  carrier  degeneracy,  and  dopant  density-dependent  dielectric 
constant  will  have  a strong  influence  on  the  mi nority-carrier  current 
flow  in  the  device  quasi -neutral  region,  and  result  in  the  modification 
of  the  nondegenerate  pn  junction  theory.  An  analysis  for  the  terminal 
behavior  and  capacitance  of  a heavily-doped  pn  diode  is  proposed  in  this 
Chapter.  The  expressions  for  the  I-V  characteristics  and  capactance  are 
modified,  which  include  the  heavy  doping  effects.  It  not  only  helps  the 
interpretion  of  the  experimental  measurements  mentioned  above,  but  also 
suggests  means  of  determining  the  bandgap  narrowing.  This  research 
effort  will  contribute  to  the  understanding  of  the  effects  of  heavy 
doping  on  the  carrier  transport  in  a heavily-doped  p-n  junction  diode. 

It  is  well  known  that  the  effect  of  bandgap  narrowing  is  to  cause 
the  bandgap  shinkage  when  dopant  density  is  increased  in  a heavily  doped 
semiconductor.  A typical  heavily-doped  pn  junction  fabricated  by  same 
material  will  thus  display  different  bandgaps  between  p-  and  n-region 
with  smaller  bandgap  found  in  the  more  heavily-doped  side.  This  bandgap 
difference  of  a heavily-doped  pn  junction  is  similar  to  that  of  a 
typical  anisotropic  heterojunction  diode.  However,  the  p-  and  n-region 
in  a heavily-doped  pn  junction  still  have  the  same  affinities.  It  gives 
us  a clue  to  study  the  effects  of  heavy  doping  by  modifying  the  theory 
of  anisotropy  of  heterojunction  diodes  [33-35]. 

In  section  4.2,  general  expressions  for  electrical  potential, 
carrier  density,  and  built-in  potential  for  a heavily-doped 
semiconductor  device  will  be  derived.  Theory  for  electrical  behavior  of 
a heavily-doped  graded  pn  junction  diode  is  then  developed.  The  results 
of  numerical  calculations  are  presented  in  section  4.3.  Section  4.4 
shows  summary  and  the  conclusions. 
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4.2  Theoretical  treatment 

4.2.1  Electrical  potential  and  carrier  densities 
In  this  section  we  cosider  a degenerate  semiconductor  with  a dopant 
concentration  that  varies  with  position.  When  the  system  is  in  thermal 
equilibrium,  the  electron  current  and  hole  current  must  separately  be 
equal  to  zero.  The  total  electron  current  consists  of  both  drift  and 
diffusion  components: 

Jn  = nnEx  + 1D„  £ <4-» 


Since  Jn  = 0,  we  may  solve  for  the  electric  field  Ex  in  terms  of  the 
electron  concentration  and  its  gradient 


E 


x 


(Dn/nn)  ^ (dn/dx) 


kT  n _1  dn 
q dn/dnn  n dx 


' - Jj1  (<vdx> 


(4.2) 


Here  we  have  used  the  generalized  Einstein  relation  for  the 
degenerate  semiconductor  derived  by  Landsberg  [13].  Since  the  Fy£ 
(Fermi-Dirac  integral  of  order  one-half)  for  the  spherical  band  in 
k -space  or  a parabolic  quantum  density  of  states  is  justified 
[5,17,36-38],  the  reduced  Fermi  energy,  nn,  for  electrons  can  be 
expressed  by  (see  Eq.  (2.1)). 
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« £n(Y  ) + F(Y  ) 
v n n7 


(4.3) 


Eq.  (4.3)  has  an  important  implication  in  that  it  allows  us  to 
derive  explicit  expressions  for  the  excess  minority  carrier  density  in  a 
heavily  doped  p-n  diode.  This  is  discussed  next. 

We  may  relate  the  electron  density  to  the  electrical  potential  via 
Eq.  (4.2).  Since  Ex  = - we  have 


at  any  point  in  an  n-type  semiconductor.  If  we  integrate  Eq.  (4.4) 
between  any  two  positions--for  example,  from  x^  to  x2,  and  through  the 
use  of  Eq.  (4.3),  we  obtain 


Since  Yn2/Ypl  = n2/nl5  we  can  rewrite  Eq.  (4.5)  in  exponential  form. 


<H  = I1  (dnn) 


(4.4) 


*n2  ' ^nl  q ^ nn2  “ 


(4.5) 


n2/n1  = exP  ^ Un2  ' * exp  CF(Ynl}  ' F(Yn2)] 


(4.6) 
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Note  that  Eq.  (4.6)  will  reduce  to  the  conventional  expression  for  the 
nondegenerate  case  if  F(Ynl)  and  F(Yn2)  are  set  equal  to  zero  (see  Fig. 
2.1).  A physical  interpretation  of  this  result  is  that  ratio  of  carrier 
densities  at  position  x2  and  x-^  is  related  not  only  to  their  potential 
difference  but  is  also  influenced  by  the  Fermi-Dirac  statistics  through 
the  function  of  F(Y). 

Similarly,  the  ratio  of  hole  densities  at  positions  x2  and  x±  may 
be  related  to  the  potential  <j>pl  and  <j>p2  by 

P2/Pl  = exP  C kT  (4>pl  ‘ *p2^  ’ exp  CF^Ypl^  " F(Yp2):i  (4'7) 

where  Yp  = p/Nv;  p is  the  hole  density,  and  Ny  is  the  effective  density 
of  the  valence  band  states. 

4.2.2  Built-in  potential 

The  expression  of  built-in  potential  for  a heavily-doped  n+p  diode 

is  derived  in  this  section.  A heavily-doped  graded  junction  n+p  diode 

with  a transition  region  modeled  to  extend  beyond  the  space-charge 
region  is  assumed  to  have  a constant  bandgap  gradient.  Variations  of 
the  dielectric  constant,  effective  mass,  carrier  lifetime,  and  mobility 
are  ignored  in  this  transition  region.  The  junction  is  said  to  be 
graded  if  the  transition  from  one  region  to  another  takes  place  over  a 
large  distance  compared  to  interatomic  distances. 

Figure  4.1(a)  shows  the  transition  region  extends  from  -Xng  to  Xpg 
and  the  space  charge  region  is  from  -Xn  to  Xp.  The  energy  band  diagram 
is  also  shown  in  Fig.  4.1(b). 
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From  Fig.  4.1,  the  built-in  potential 


can  be  expressed  by 


G _ * * 

♦bi  = *nl  ‘ V 

6E 

= +nl  + < | [Eg  - <xpg  - xp>]  - +p2l  (4-8> 

where  L = xpg  + xng,  6Eg  = - E^J;  E^  and  E^J  are  the  bandgap  of  p-  and 

n-region,  respectively.  From  Eqs.  (2.10)  and  (2.11),  we  obtain 


G 1 |-rP  ^g  , n , * * 

bi  = q CEg  L xpg  - xp)]  + <+nl  “ *p2} 


■rC£n  (n^)3+^CF(Yn)+F(Yp)]+IE;G 


(4.9) 


and 


EgG  =?[(Eg  - 4EgP)  'T2  (xpg  - xp)] 


(4.10) 


where  F(Yn)  and  F(Yp)  are  given  by  Eq.  (2.3)  and  may  be  considered  as 
corrections  in  the  Fermi-Dirac  statistics  for  the  degenerate  region.  As 
shown  in  Fig.  2.1,  both  the  values  of  F(Yn)  and  F(Yp)  will  approach  to 
2 x 10"3  and  1 x 10'3  when  doping  density  is  equal  to  1017  cm'3.  Note 
that  the  second  term  in  the  right-hand  side  of  Eq.  (4.10)  is  formed  from 
the  graded  reduction  of  bandgap  in  the  transition  region  of  the 
junction . 
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4.2.3  Graded  junction  diode 

Since  a heavily  doped  p-n  junction  diode  is  treated  as  a device 
with  a linearly  varying  bandgap  in  the  transition  region,  additional 
terms  in  the  transport  equation  in  addition  to  that  of  the  conventional 
drift  and  diffusion  will  arise  from  the  variation  of  the  energy  band 
edges  and  from  the  spatial  dependence  of  the  electron  and  hole  density 
of  states,  where  the  rigid  band  approximation  is  assumed.  We  have  also 
made  the  fixed-field  approximation  which  is  valid  for  low-level 
injection  and  depletion  approximation  to  facilitate  the  solution  of 
Poisson's  equation.  The  model  also  assumes  uniform  doping  in  both  n and 
p regions.  The  electron  and  hole  current  densities,  in  the  rigid  band 
case,  may  then  be  written  as 


electron  and  hole  diffusion  coefficients,  respectively;  <f>  is  the 
electrostatic  potential,  and  Eg  denotes  the  bandgap.  Since  Jp  = 0 for 
the  quasi -neutral  p-region  in  thermal  equilibrium,  we  have,  for  low- 
level  injection 


J = -ny  (q<|>)  + qD 
n n 3x  M n 3x 


(4.11) 


(4.12) 


Here,  Up  and  yn  are  the  hole  and  electron  mobilities;  DR  and  Dp  are  the 


(4.13) 
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For  uniform  doping,  the  injected  electron  and  hole  currents  may  thus  be 
expressed  by 


i ix  3n  . 
Jn  = kTl,n  W + 


dE 

nu  x— 2- 
n dx 


n 3n  , n 1 
= ^Dn  97  + ‘'"“n  r 


X > X 


» dE 

Jp  = ‘ kTyp  ‘ pyp  d^ 


- qD  •J—  - qpD  -j — , x ^ x 
H p 3x  p Lg  * n 


(4.14) 


(4.15) 


where  Lg  is  assumed  to  be  constant  in  the  transition  region.  The 
continuity  equations  for  the  injected  electrons  and  holes,  under  steady 
state  condition,  can  be  written  as 


1 

q 


^n 

dx 


n 

T 

n 


1 

q 


(4.16) 


(4.17) 


where,  xn  and  Tp  are  the  minority  electron  and  hole  lifetimes, 
respectively.  Here  xn  and  xp  are  no  longer  only  depending  on  the 

Schockl ey-Hall -Read  (SHR)  recombination  through  recombination  centers  in 
the  degenerate  semiconductor,  but  trap-assisted  Auger  and  band-to-band 
Auger  recombination  are  important  recombination  processes  which  should 
be  included.  p°  (or  n°)  denotes  the  thermal  equilibrium  carrier 

A 

n and  p are  the  excess  electron  and  hole  densities. 


density. 
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respectively.  For  the  boundary  conditions  that  An  = Ap  = 0 at  x = + “, 
and  n,  p and  their  gradients  are  continuous  everywhere,  the  current  and 
continuity  equations  for  the  minority  carriers  are  solved  by  [33]. 
Thus,  the  electron  and  hole  densities  are  given  respectively  by 


n(x)  = { 


+ c^2'  + exp  [-(x  - xpg)/Lg]  , xp  < x < xpg 

-x/Ln  (n?/ 

C3e  + Na  ’ x > xpg 


(4.18) 


, N % 2 

«iX  apx  (n.  ) 

B,e  1 + B0e  + * exp  [-(x  + x )/L  ],  -x  < x < -x 


1 


P(x)  = { 


'D 


x/L  (nje)2 

B3e  + -N^->  x < ~xng 


(4.19) 


where  L = (D  t )1//2  L = (D  t ) 1//2  are  the  minority  carrier  diffusion 
n n n p p p 

length;  n^e  and  n?e  are  defined  by  Eq.  (2.27)  and  (2.28).  Other 
constants  are  defined  by 


{1  + [1  + 

(2Lg/Ln)2]1/2} 

(4.20a) 

{1  - [1  + 

(2Lg/Ln)2]1/2l 

(4.20b) 

C,  = n(x  )ek2Xpge  ^Sk,,  + f-)/D 
1 p n 


(4.20c) 
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and 


„ k.x  -x  /L  , 

-n(x)e  1 P9e  ^ "(k,  ♦ f-)/D 

p n 

( 4 ,20d ) 

- k1Xr,^  k rtX 

n(Xp)e  1 P9e  2 P9(k?  - k 1)/D 

(4.20e) 

-x  /L  k,x  k0x  , 

e P9  " . [e  1 P • e 2 p9(k  + f-) 

Ln 

k 0x  k , x . 

- e 2 Pe  1 P9(k  + 1 )] 

n 

( 4 . 20f ) 

- or-ll  + [1  + (2L q/L  )]1/2J 
g y p 

(4.21a) 

- kt  (i  - [i  - <2yLp)]1/2} 

(4.21b) 

P<-x„)e  2n9e  "9  <>(a,  - f-)/G 

P 

(4.21c) 

~ t \ ~alXn9  -X ng^ *~p / 1 %/r 

-p(-x  )e  ye  y H(ax  - — )/G 

(4.21d) 

-a, x_  -a0x 


A '•*  1 A VArtA 

p ( -x n ) e 9e  ng(a2  - c*1)/G 


(4.21e) 


e "9  P[e  2 " 


-a,x  -oux 


-X  /L  -a„x„  ~a1X„„  i -vaoA  , 

• e 1 "9(f-  0l)  - e 1 n • e 2 "9(-L  - «,)] 
P P 


( 4 . 21 f ) 


By  substituting  Eqs.  (4.18)  and  (4.19)  into  Eqs.  (4.14)  and  (4.15),  we 
obtain  the  electron  and  hole  currents  at  the  space  charge  region 


boundaries  as 
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j„(xp)  = -dV(xP>xi 


(4.22) 


Jp(-Xn>  = -q°pP<-xn>J-2 


(4.23) 


where 


X - = - e 


klxp  C1 


n (x  ) 
P 


k 0x  C 


(kx  + ±-)  - e 2 p (k2  + j-) 


n(x  ) 
P 


(4.24) 


and 


-a-|  X , -CUX  . 

x2  - e (ai  + T- ) + e (a?  + r— ) 

P(-xn)  1 9 p(-xn)  2 Lg 


(4.25) 


By  assuming  no  recombination  in  the  space-charge  region,  the  total 
forward  current  density  is  given  by 


Jtot  = c'V('Xn)12  + <'Dnntxp,Al 


(4.26) 


The  minority  carrier  densities  at  the  space-charge  region 
boundaries  may  be  obtained  by  using  law  of  the  junction  as 


n (x  ) = n°(x  ) exp  (qV  /k T) 
P p a 

(nfL(x™))2 


— exp  [(x  - x )/L  ] • exp  (qV  /kT) 


pg  p g- 


(4.27) 
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p(-xn)  = P°(-xn)  exp  (qVa/kT) 


• exp  [-  (xpg  - xn)/L  ] • exp  (qVa/kT) 


(4.28) 


Thus,  the  diode  current  density  in  Eq.  (4.26)  can  be  expressed  as 


It  is  noted  that  Eq.  (4.29)  will  reduce  to  the  conventional  diode 
current  density  expression  in  the  limit  of  1/L  = 1/Lg  = 0.  The 

additional  exponential  terms  appeared  in  Eq.  (4.29)  are  due  to  the 
effect  of  bandgap  narrowing. 

The  space  charge  region  boundaries,  Xn  and  Xg  may  be  obtained  as 
functions  of  the  applied  voltage  by  solving  Poisson's  equation  as 
described  in  section  2.3.3.  We  obtain  the  voltage  drop  across  the 

space-charge  region  as 


( n . ) 

e«p[(xn  - xng)/Lg]  + «ID„X1 — ^ — exp|;(xpg  - Xp)/Lg]) 


x [exp  (qV  /k T ) - 1] 

a 


= J • [exp(qV  /k T)  - 1] 

U a 


(4.29) 


(4.30) 
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Since  the  total  charge  per  unit  area  in  the  space-charge  region  is  Q^qR 
= qNDxp  - qNAXp,  and  I QsCR I <<:  c'^Dxn  or  qNAxp5  we  may  assume 


Vn  “ Kflxp 


(4.31) 


The  space  charge  boundary,  Xp  and  Xp  for  a heavily-doped  pn  graded 
junction  can  be  obtained  by  using  Eq.  (4.30)  and  (4.31) 


o i - v 1/2 

r (1_  + Ik] a) 

qNA  L(c2  C > 


.6 


x_  = 


2? 


n qN 


[( 


1 , ♦bi  - V/2 


D C 


V 


& 


(4.32) 


(4.33) 


where 


z,  = (qe  e NnNn)/2(e  N,  + e Nn) 
VM  sp  sn  A D'  ' sp  A sn  D' 


(4.34) 


r = 6Vq2NAL 

g 


(4.35) 


Thus,  the  depletion  layer  width  can  be  written  as 


,2? 


1 . 1 . «b1  ' V/2  1 


xd  ’ ( q)(NA  + Nd>[<c?  + 


-] 


(4.36) 


where  <j>^  is  the  built-in  potential  given  by  Eq.  (4.9). 

The  transition  region  boundaries,  Xpg  and  xpg,  are  experimentally 
determined  by  the  plot  of  vs . Vfl.  It  may  be  depicted  as  follows: 

c 

The  capacitance  per  unit  area  is  defined  by  C = |dQ/dVa|,  where  Q 


represents  the  charge  per  unit  area  in  half  of  the  junction  region.  We 


can  write 
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C 


dx 

qNA  dV2, 


(4.37) 


By  substituting  Eq.  (4.30)  and  (4.31)  into  Eq.  (4.37),  we  obtain 


C = 


(xn/Esn>  + <xp/esp>  + r 


(4.38) 


Substituting  Eq.  (4.31)  and  (4.32)  into  (4.38),  we  have 


1 = ( 1 + ^ " s 


(4.39) 


1_ 

where  ^2  is  due  to  the  bandgap  narrowing  effect  imposed  on  the  electric 
9 

field  at  the  edges  of  the  space  charge  region.  From  Eq.  (4.39),  a plot 
of  Kj-  vs.  Vg  gives  a straight  line  with  an  extrapolated  interception 

c 

with  voltage  axis  at 


+ ?/Cg  (4-«» 

Thus,  Xpg  can  be  determined  from  Eq.  (4.40)  and  (4.9)  by  using  the 
dopant-density  dependent  bandgap  narrowing  data  and  assuming  Xpg  = xpg 

” 2* 


4.3  Results  and  discussion 

Calculations  of  the  effects  of  heavy  doping  on  the  built-in 
potential  and  electron  current  density  for  a heavily-doped  Si  n+p  graded 
junction  were  carried  out  by  assuming  a constant  bandgap  gradient  and 
using  the  expressions  described  in  section  4.2.  Table  4.1  shows  the 
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calculated  built-in  potential  of  a heavily-doped  Si  n+p  graded  junction 
with  and  without  including  the  heavy  doping  effects,  and  assuming  xng  = 
Xpg  = L/2  for  different  transition  region  widths.  The  doping  density  of 
n+  side  is  assumed  to  be  1 x 10^  cnT^,  and  the  acceptor  density  on  the 

1 o o 

lightly  doped  side  is  5 x 10  cm”  . The  bandgap  narrowing  data  used  in 
these  calculations  are  reported  by  Neugroschel  et  al . [7].  From  Table 

p 

4.1,  it  is  shown  that  the  modified  built-in  potential,  <()^.,  calculated 
from  Eq.  (4.9)  assuming  L = 0.05,  0.1,  0.2,  and  0.4  ym  and  without 

including  the  effect  of  dopant  density  dependent  dielectric  constant, 
are  about  110  to  140  mV  smaller  than  those  calculated  from  conventional 

g 

expression.  Values  of  4^  calculated  from  Eq.  (4.9)  including  the 
effect  of  dopant  density  dependent  dielectric  constant  (assuming  e$p  = 
32.9)  are  also  listed  in  Table  4.1.  The  values  are  about  90  to  140  mV 
smaller  than  those  calculated  from  conventional  expression.  Thus,  the 
heavy  doping  effects  such  as  carrier  degeneracy  and  bandgap  narrowing 
have  an  important  influence  on  the  value  of  built-in  potential  of  a 

heavily-doped  Si  pn  graded  junction.  Note  that  the  effect  of  dopant- 
density  dependent  dielectric  constant  will  signficantly  increase  the 

Q 

values  of  ^ in  the  narrow  transition  region.  As  for  a wider 
transition  region,  the  effect  of  dopant  density  dependent  dielectric 
constant  can  almost  be  ignored. 

The  results  of  the  calculated  electron  current  density  for  a 

heavily-doped  Si  n+p  graded  junction  are  shown  in  Figs.  4.2  and  4.3,  in 

which  the  doping  densities  of  the  n+  and  p side  are  assumed  to  be  1 x 
10^0  and  5 x 10^  cnT^,  respectively.  The  corresponding  minority 

carrier  diffusivities , Dp  and  DR,  are  assumed  to  be  1.0  and  5.0  cm^/sec. 
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respectively;  where  the  minority  carrier  diffusion  lengths,  Lp  and  Ln, 
are  equal  to  0.25  and  21  ym  [32,21],  respectively.  Fig.  4.2  shows  the 
electron  current  density  vs.  applied  forward  voltage  of  a heavily-doped 
Si  n+p  graded  junction  with  a constant  bandgap  gradient  and  e = 

b p 

11.7.  Solid  line  is  the  values  calculated  from  Eq.  (4.22)  with  L = 0.05 
ym,  dash  line  with  L = 0.1  ym,  and  dash-dot  line  with  L = 0.4  ym, 
respectively;  where  xpg  = xpg  = L/2  is  assumed.  The  results  show  that  a 
smaller  transition  region  width  will  result  in  smaller  electron  current 
density  in  the  quasi-neutral  n-region.  Fig.  4.3  also  shows  the  electron 
current  density  vs.  applied  forward  voltage  for  the  same  heavily-doped 
Si  n+p  graded  junction  with  and  without  including  the  heavily-doped 
effects  for  L = 0.1  ym.  Solid  line  is  the  values  calculated  from  Eq. 
(4.22)  with  e$p  = 11.7,  dash-dot  line  with  £$p  = 32.9,  and  dash  line  is 
obtained  from  the  conventional  expression,  respectively.  The  results 
show  that  the  effects  of  heavy  doping,  especially  the  bandgap  narrowing 
and  the  carrier  degeneracy,  have  a significant  influence  on  the  electron 
current  density  at  the  quasi -neutral  region  of  a heavily-doped  Si  pn 
graded  junction.  This  is  probably  because  the  bandgap  gradient  affects 
the  transport  of  the  injected  minority  carriers  by  influencing  the 
carrier  distribution  in  the  quasi -neutral  regions.  The  current  density, 
therefore,  depends  on  the  width  of  the  transition  region  outside  the 
space  charge  region.  Since  this  width  is  determined  by  the  width  of  the 
space-charge  region  on  either  side,  J will  also  depend  on  the  applied 
voltage  through  xn  and  xp.  and  X2  are  the  only  weak  functions  of 
applied  voltage,  while  exp(xn/Lg)  is  the  dominant  factor. 
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The  dependence  of  current  density  on  bandgap  gradient  is  due 
primarily  to  the  fact  that  thermal  equilibrium  pn-product  has  to  be 
modified  for  a semiconductor  with  position-dependent  band  structure. 
Thus,  the  dark  saturation  current  in  a pn  junction  as  well  as  the  common 
emitter  gain  of  a bipolar  transistor  are  directly  affected  by  the 
modified  pn-product. 

4.4  Summary  and  conclusions 

A theoretical  approach  of  studying  the  heavy  doping  effects  on  the 
built-in  potential  and  current  density  in  a heavily-doped  Si  n+p  graded 
junction  has  been  developed  in  this  chapter.  A closed  form  expresson 
for  the  I-V  characteristics  is  presented  by  modifying  the  theory  of 
anisotropy  heterojunction  diode  [33-35],  in  which  the  depletion 
approximation  and  the  rigid  band  model  are  assumed.  The  expression  is 
valid  for  uniform  doping  and  low-level  injection  cases.  A transition 
region  is  modeled  to  extend  beyond  the  space  charge  region  with  a 
constant  bandgap  gradient.  Numerical  calculations  of  the  heavy  doping 
effects  such  as  carrier  degeneracy,  dopant  density  dependent  dielectric 
constant,  and  bandgap  narrowing  were  carried  out  by  employing  the  models 
described  in  Section  2.2.  The  effects  of  heavy  doping  and  variations  of 
transition  region  width  on  the  built-in  potential  and  current  density  of 
a heavily-doped  Si  n+p  graded  junction  are  theoretically  analyzed  in 
this  chapter. 
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TABLE  4.1 


A comparision  of  the  calculated  built-in  otential  of  a heavily-doped  Si 
pn  graded  junction  with  and  without  including  the  heavy  doping  effects 
for  different  transition  regin  width. 


Built-in  Potential 
(Calculated  at 
Va  = 0.1  Volt) 

L 

(pm) 

0.05 

0.1 

0.2 

0.4 

^bi 

(Gsp 

= 11.7) 

0.989 

0.971 

0.963 

0.958 

^bi 

( esp 

= 32.9) 

1.01 

0.984 

0.969 

0.962 

( Nondegenerate 

1.10 

expression) 
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(b) 


Figure  4.1  (a)  Model  of  the  N - P heavily-doped  graded  junction; 

Vj  is  the  voltage  across  the  space  charge  region 
and  L is  the  width  of  the  transition  region. 

(b)  The  corresponding  energy  band  diagram. 


Electron  Current  Density  Jn  ( A/cm2) 
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Figure  4.2  Electron  current  density  vs.  applied  forward  bias  with 
transition  region  width  as  a parameter  for  a heavily- 
doped  pn  graded  junction.  Solid  line  is  the  values 
calculated  from  Eq.  (4.22)  with  L = 0.05  jjm,  dash  line 
with  L = 0.1  jjm,  and  dash-dot  line  with  L = 0.4  pm,  res- 
pectively; where  xng  = xpg  = L/2  is  assumed. 


Electron  Current  Density  Jn(A/cm2) 
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Applied  Voltage  Va  (Volt) 


Figure  4.3  Electron  current  density  vs.  applied  forward  bias  for 
a heavily-doped  Si  pn  graded  junction.  Solid  line  is 
the  values  calculated  from  Eq.  (4.22)  with  dielectric 
constant  equal  to  11.7,  dash-dot  line  with  dielectric 
constant  equal  to  32.9,  and  dash  line  is  obtained  from 
the  conventional  expression,  respectively. 


CHAPTER  V 

EFFECTS  OF  HEAVY  DOPING  ON  THE  BUILT-IN  POTENTIAL 
DEPLETION  LAYER  WIDTH,  DEPLETION  CAPACITANCE,  AND 
I-V  CHARACTERISTICS  OF  A HEAVILY-DOPED  HETEROJUNCTION  DIODE 

5.1  Introduction 

Recently,  two  new  epitaxial  technologies,  namely,  the  molecular 
beam  epitaxy  (MBE)  and  the  metal -organic  chemical  vapor  deposition 
(MOCVD),  have  emerged  into  semiconductor  technologies.  They  offer  the 
promise  of  making  highly  advanced  heterostructures  routinely 

available.  Among  the  different  combinations  of  material  being  used  to 
produce  a heterojunction,  the  GaAlAs-GaAs  system  appears  to  be  the  most 
promising  candidate.  This  is  due  to  the  fact  that  Ga ( i_x ) A1 xAs  lattice 
closely  matched  to  that  of  GaAs,  and  the  heterojunction  structure  can  be 
formed  without  severe  interface  state  effects. 

While  many  kinds  of  devices  will  benefit  from  these  emerging 
technologies,  the  bipolar  transistors  and  high  electron  mobility 
transistors  (HEMT)  (or  Modulation-doped  field  effect  transistors 
(MODFET))  using  III-V  compound  semiconductor  materials  will  benefit  most 
from  the  MOCVD  and  MBE  technologies. 

The  underlying  physics  is  the  use  of  energy  gap  variations  besides 
electric  field  to  control  the  forces  acting  on  electrons  and  holes,  it 
permits  a re-optimization  of  doping  levels  and  geometries,  leading  to 
high  speed  devices.  The  potential  advantages  of  heterojunction  bipolar 
transistors  (HBJT)  have  been  described  [39],  For  example,  a high 
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injection  efficiency,  a decrease  in  base  resistance,  a reduction  of  the 
emitter  current  crowding,  and  an  improvement  in  frequency  response  (high 
gain  bandwidth  product)  have  become  very  promising  possibilities  in  the 
near  future.  In  addition,  the  doubl e-heterostructure  transistor  with 
both  wide-gap  emitters  and  collectors  offer  additional  advantages.  It 
exhibits  a better  performance  under  saturated  operation,  and  their 
emitter  and  collector  may  be  interchanged  by  simply  changing  biasing 
conditions,  greatly  simplifying  the  architecture  of  bipolar  IC's. 

HEMT  have  been  intensively  studied  by  many  researchers  at  different 
laboratories.  This  field-effect  transistor  is  fabricated  by  a 
heterojunction  with  a doped  AlGaAs  layer  to  an  adjacent  undoped  GaAs  in 
such  a way  that  the  mobilities  and  velocities  of  the  charge  carriers 
moving  in  the  undoped  gallium  arsenide  can  be  modulated  and  optimized. 
The  device  also  accepts  and  releases  considerable  amounts  of  charges 
rapidly,  and  it  dissipates  very  little  energy  in  switching  operation. 
Thus  it  is  also  expected  to  be  the  very  promising  candidate  for  the 
needs  of  supercomputers  [40], 

Both  the  HBJT  and  HEMT  contain  heavily-doped  heterojunction,  thus 
it  becomes  very  important  to  study  the  effects  of  heavy  doping  on  the 
built-in  potentail , depletion  layer  width,  depletion  capacitance,  and 
I-V  characteristics  in  a heavily-doped  heterojunction  diode.  Since  the 
effective  density  of  conduction  band  states  of  GaAs  and  A1 Q 3GaQ  yAs  are 
about  4.7  x 10^  and  7.0  x 10^  respectively  at  T = 300  K,  carrier 
degeneracy  will  be  the  dominant  heavy  doping  effect  in  the  heavily-doped 
GaAs-Al GaAs  heterojunction  diode.  While  the  effect  of  bandgap  narrowing 
will  be  less  important  due  to  the  very  small  electron  effective  mass  in 
GaAs  and  AlGaAs. 
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Both  MBE  and  MOCVD  technologies  are  capable  of  growing  epitaxial 
layers  with  high  crystalline  perfection  and  purity.  Highly  controlled 
doping  levels  up  to  1019  impurities  per  cm3  and  higher  can  be 
achieved.  In  addition,  highly  controlled  changes  in  doping  level  are 
possible  during  growth  without  interrupting  the  epitaxial  growth,  and 
with,  at  most,  a minor  adjustment  in  growth  parameters  [39].  The  doping 
may  be  changed  either  gradually  or  abruptly  due  to  the  comparatively  low 
growth  in  temperatures.  Effects  of  heavy  doping  on  a heavily-doped 
graded  and  abrupt  heterojunction  will  thus  be  depicted  in  the  following 
sections.  In  section  5.2,  the  expression  for  built-in  potential  is 
derived,  and  expressions  for  the  depletion  layer  width  and  depletion 
capacitance  are  reformulated.  Results  of  numerical  calculations  are 
described  in  section  5.3.  Section  5.4  presents  a summary  and 
concl usions . 

5.2  Graded  heterojunction 

5.2.1  Built-in  potential 

In  addition  to  the  assumptions  made  for  a heavily-doped  Si  pn 
junction  diode  as  described  in  Chapter  IV,  we  also  assume  that  the 
change  in  band  structure  due  to  strain  is  negligible. 

The  energy  band  diagram  before  formation  of  a pn  GaAs-Al  Ga-i  As 
heterojunction  diode  is  shown  in  Fig.  5.1.  Where  x and  $ denote  the 
electron  affinity  and  work  function  of  semiconductor,  respectively. 
Fig.  5.2  shows  the  schematic  representation  of  a GaAs-Al xGa1_xAs  pn 
graded  junction  in  thermal  equilibrium.  We  may  obtain  the  built-in 
potential  across  the  space  charge  region 
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where  AX  is  the  electron  affinity  difference  across  the  space  charge 
boundaries.  A hyperbolic  tangent  is  usually  a good  approximation  for  X 
[41],  i.e., 

X (x ) = (6X/2)  tanh  (mx/L)  + b (5.2) 


Here,  m and  b are  constants,  SX  = Xj  - X2,  and  L = xpg  + xpg.  If  the 
width  of  the  space  charge  region  is  small  compared  to  that  of  the 
transition  region,  X(x)  can  be  approximated  by  its  Taylor  series 
expansion  up  to  the  linear  term.  Thus, 


and 
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(5.3) 

(5.4) 


For  a heavily-doped  graded  heterojunction , we  may  use  Eq.  (2.1),  (2.2) 
and  (5.4),  and  modify  Eq.  (5.1) 
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kT  r„_  AV-,  , kT 


q - [*"  (NX)]  + q CF(Yn}  + F(Yp)]  + EgG^  (5*5a) 


V c 


Where,  F(Yp)  and  F(Y  ) are  defined  by  Eq.  (2.3),  and 


6E 


EgG  = [Egl  - 4Egl  + T3'  (xpg  ' xp>  ' <xn  + xp>  f ™]  (5-5b> 


here,  AEg-^  is  the  bandgap  narrowing  in  the  heavily-doped  p region. 


5.2.2  Depletion  layer  width  and  depletion  capacitance 
The  depletion  layer  boundaries,  xn  and  Xp,  can  be  obtained  by 
solving  Poisson's  equation 


P - ND  + Na) 


d<(> 
" dx 


d£ne 
s 

dx 


(5.6) 


Since  accurate  numerical  solution  for  an  n+-p  structure  with  a Gaussian 
doping  profile  indicates  that  even  a substantial  variation  of  e has 
little  effect  on  the  solution  for  <f> ( x ) [33],  Eq.  (5.6)  can  be  simplified 
by  assuming  es  to  be  constant  in  the  region  of  interest  and  by  making 
the  depletion  approximation.  The  boundary  conditions  for  the  quasi- 
neutral region  are  given  by 
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Solving  the  Poissons  equation  under  these  conditions  we  obtain  the 
voltage  drop  across  the  space-charge  region 
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(5.9) 


Solving  Eq.  (5.9)  and  (4.31)  simultaneously  and  using  Eq.  (5.3),  we 
obtai n 


where 
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(5.11) 
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(5.12b) 


Thus  the  depletion  layer  width  for  a heavily-doped  graded  hetero junction 
can  be  expressed  by 


‘dG  = <r>%  + ^<c2 
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where,  <()bi  is  given  by  Eq.  (5.5). 
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The  depletion  capacitance  can  also  be  written  as 


C 


dx 

qNA  dV1 
a 


(5.14) 


O 

where,  the  term  1/Cg^  is  due  to  the  changes  of  bandgap  and  electron 

affinity  at  the  edges  of  the  space  charge  region.  Using  the  dopant- 

density  dependent  bandgap  narrowing  data  and  assuming  x = x = L/2, 

r y y 

the  transition  region  boundaries,  xng  and  Xpg,  can  also  be  determined  by 
the  plot  of  1/C2  vs.  Vg  and  Eq.  (5.5). 


5.3  Abrupt  heterojunction 

An  abrupt  heterojunction  is  distinguished  from  a graded 
heterojunction  according  to  the  width  of  the  transition  region  near  the 
interface  of  the  two  semiconductors.  For  a heterojunction  to  be  abrupt, 
the  transition  from  one  material  to  the  other  is  completed  within  a few 
atomic  distances.  For  this  case,  the  material  parameters  are  assumed  to 
change  abruptly  at  the  interface.  We  assume  this  interface  is  formed 
between  semiconductors  with  a good  lattice  match,  e.g.,  GaAs-A^Ga-^As, 
so  that  the  effects  of  lattice  strain,  interface  states,  and  dipoles 
[34]  can  be  neglected  (i.e.,  no  interfacial  charge  at  the  interface). 

5.3.1  Built-in  potential 

Fig.  5.3  shows  the  schematic  representation  of  GaAs-Al xGa^_xAS 
heterojunction  in  thermal  equilibrium.  Here,  ^ is  the  total  built-in 
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potential,  <f>Qp  is  the  built-in  potential  on  the  P-side,  and  on  the 
N-side  is  4>DN . The  total  built-in  potential  across  the  space  charge 
region  is 


*D  = *n  ' % 
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Eq.  (5.15)  can  be  written  as 
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Conventionally,  the  built-in  potential  for  an  abrupt  nondegenerate 
doping  heterojunction  is  given  by  [33] 
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For  a heavily-doped  abrupt  pn  heterojunction,  Eq.  (5.17)  can  be  modified 
using  Eq.  (2.1)  and  (2.2) 
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Yn  = "/Nc2  * Yp  = P/Nvl  * 

EgA  = 2 (Egl  + Eg2^  " \ (AEgl  + AEg2) 


(5.19a) 


(5.19b) 


here,  AEgl  and  AEg2  are  the  bandgap  narrowing  due  to  the  effect  of  heavy 
doping. 


5.3.2  Depletion  layer  width  and  depletion  capacitance 
Since  the  electric  flux  density  D = e$E  is  continuous  at  the 
junction  interface  (i.e.,  e$nE^  = espEp  at  X = 0)  and  <f>Q  = <J>q^  + <j>gp,  we 
can  obtain  the  depletion  region  boundaries,  xn  and  xg,  for  a heavily- 
doped  abrupt  heterojunction  by  using  the  nondegenerate  formulations  from 
[42]  and  replacing  the  conventional  built-in  potential  with  Eq.  (5.19) 
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The  depletion  capacitance  thus  is  given  by 


(5.23) 


By  plotting  1/C2  as  a function  of  Va,  Eq.  (5.23)  shows  that  the 

O 

intercept  of  1/C  = 0 is  4>q . Now  by  measurement  of  capacitance  per  unit 


area,  AEg  (=  AEgl  + AEgg)  may  be  found  from  Eq.  (5.19)  if  the  AEC 
(=  - X2),  Fermi  levels,  dopant  densities,  and  intrinsic  energy  gaps 

on  each  side  of  the  hetero junction  are  known.  Due  to  the  very  small 
electron  density  of  state  effective  mass,  the  calculated  values  of  AEg 
will  be  mainly  contributed  by  AEgl.  Thus,  AEgl  may  be  separately 
determined  from  such  calculations. 

5.3.3  I-V  characteristics 

In  the  hetero junction  case,  the  diffusion  current  is  dominated  by 
injection  from  the  wider  energy  gap  semiconductor  into  the  narrow  energy 
gap  semiconductor.  There  are  three  assumptions  for  the  ideal  current- 
voltage  model:  1)  the  space  charge  region  depletion  layer  has  abrupt 

boundaries  so  that  outside  the  depletion  region  the  semiconductor  is 
assumed  to  be  neutral;  2)  the  injected  minority  carrier  densities  are 
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assumed  small  compared  to  the  majority  carrier  densities;  3)  there  is  no 
generation  and  recombination  current  in  the  depletion  region  so  that  the 
electron  and  hole  currents  are  constant  through  the  depletion  region. 

Thus,  from  Eq.  (4.6)  and  Fig.  5.3,  we  may  obtain  the  electron 
current  flow  at  X = -xp 
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dn 
n dx 
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(5.24) 


This  flow  of  electrons  in  the  negative  x direction  corresponds  to  a 
positive  current  in  the  positive  x direction.  Similarly,  for  holes  on 
the  N-side,  the  hole  current  density  is  given  by: 


qDDNA  qV  /kT 

Jp  = Tp — exp  + 6Ev)/kT)  + p(Yp)^  * (e  - 1) 

= Jpo  * NA  * exp  + 6Ev)/kT)  + p(Yp)3  (5.25) 


The  influence  of  the  heavy  doping  effects  on  a heavily  doped 
heterojunction  transisitor  can  be  illustrated  by  taking  the  ratio  Jn/Jp: 
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Note  that  for  the  example  of  Fig.  5.3  with  X = 0.3,  where  6EC  + 6EV  = 
(1.798  - 1.424)  = 0.374  eV,  the  exponential  term  in  Eq.  (5.26) 
dominates,  and  |Jn|  >>  |J  |.  Therefore  the  diffusion  current  is  due  to 
injection  of  majority  carriers  from  the  wider  energy  gap  side  into  the 
narrower  energy  gap  side.  The  total  current  density  may  then  be 
generally  approximated  by  Jt  = Jp. 

5.4  Results  and  discussion 

Calculations  of  the  effects  of  heavy  doping  on  the  total  built-in 
potential,  depletion  layer  width,  and  hole  current  density  for  a 
heavily-doped  pn  abrupt  heterojunction  were  carried  out  by  using  the 
expressions  derived  in  Section  5.3.  Fig.  5.4  and  5.5  show  the  values  of 
function  F(Y)  vs.  hole  or  electron  density  in  the  heavily-doped  GaAs  and 
Alo.3Gag.7As  at  T = 300  K,  respectively,  where  F(Y)  is  calculated  from 
Eq.  (2.3).  Solid  lines  are  values  for  electrons,  dash  lines  for 
holes.  Values  of  F(Y)  vary  from  8 x 10"^  to  6.13  and  5.2  x 10"^  to  4.8 
for  n-type  GaAs  and  Alg  gGag  yAs  respectively  for  dopant  densities  in 
the  range  of  1 x 10*7  to  1 x 10^  cm-^. 

The  numerical  results  for  the  total  built-in  potential  of  a 
heavily-doped  np  abrupt  heterojunction  are  shown  in  Figs.  5.6  and  5.7. 
Fig.  5.6  shows  the  total  built-in  potential  vs.  emitter  dopant  density, 
Ng,  for  a heavily-doped  abrupt  Alg  gGag  yAs-GaAs  emitter-base 
heterojunction  of  the  wide-bandgap  npn  transistor  at  T = 300  K,  in  which 
the  base  dopant  density,  N^,  is  assumed  to  be  1 x 10^  cm"^.  Solid  line 
is  the  predicted  values  by  using  Eq.  (5.19),  dash  line  is  calculated  by 
using  the  conventional  theory  from  Eq.  (5.18).  Fig.  5.7  shows  the  total 
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built-in  potential  vs.  base  dopant  density,  NA,  for  the  same  heavily- 
doped  abrupt  A1  g^Gag^As-GaAs  heterojunction  at  T = 300  K,  and  the 
emitter  dopant  density,  Ng,  is  assumed  to  be  1 x 10*7  cm-3.  Solid  line 
is  the  calculated  values  by  using  Eq.  (5.19),  dash  line  is  calculated  by 
using  the  conventional  theory  from  Eq.  (5.18).  Thus,  the  heavy  doping 
effects  of  carrier  degeneracy  and  bandgap  narrowing  will  also  have  a 
significant  influence  on  the  total  built-in  potential  of  a heavily-doped 
pn  abrupt  heterojunction.  It  is  noted  that  owing  to  the  effective 
masses  of  electron  in  GaAs  and  GaAlAs  are  very  small  , the  effect  of 
bandgap  narrowing  in  n-type  GaAs  and  AlGaAs  is  negligible.  On  the  other 
hand,  carrier  degeneracy  is  the  dominant  heavy  doping  effect  in  the 
heavy-doped  GaAs-AlGaAs  hetero junction  devices. 

Fig.  5.8  shows  the  depletion  layer  width  vs.  emitter  dopant 

density,  Ng,  for  a heavily-doped  abrupt  Alg  3Gag  yAs-GaAs  heterojunction 

at  T - 300  K,  in  which  base  dopant  density,  NA,  is  assumed  to  be  1 x 
20  -3 

10  cm"  . Solid  line  is  the  calculated  values  by  using  Eq.  (5.22), 

dash  line  is  calculated  by  using  the  conventional  nondegenerate 
expression  as  shown  in  [42].  The  results  show  that  the  effects  of 
carrier  degeneracy  and  bandgap  narrowing  have  little  influence  on  the 
depletion  layer  width. 

The  hole  current  density  vs.  emitter  dopant  density  at  T = 300  K is 
shown  in  Fig.  5.9;  where  the  base  dopant  density,  NA,  is  assumed  to  be  1 
x 1017  cm"3.  Solid  line  is  the  calculated  values  using  Eq.  (5.25),  dash 
line  is  calculated  by  using  the  conventional  theory  for  abrupt 
heterojunction  diodes  [42],  Note  that  the  value  of  conduction  band 
discontinuity  used  in  these  calculations  is  given  by  [42] 
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6Ec  = (0.85  ± 0.03)  • (E^  - E^) . (5.27) 

From  Fig.  5.9,  it  is  clearly  shown  that  carrier  degeneracy  and  bandgap 
narrowing  have  an  important  effect  on  the  current  density.  The  results 
show  significant  difference  from  those  predicted  by  the  nondegenerate 
expressions.  The  effects  of  carrier  degeneracy  and  bandgap  narrowing  on 
heavily-doped  heterojunction  transistor  current  gain  are  shown  in  Fig. 
5.10,  where  the  emitter  dopant  density,  ND,  is  assumed  to  be  1 x 1017 
cm  . Solid  line  is  the  calculated  values  by  using  Eq.  (5.26);  dash 
line  is  calculated  by  using  the  conventional  theory  as  shown  in  [42]. 
Thus,  it  is  clearly  shown  that  carrier  degeneracy  have  significant 
influences  on  the  current  gain  of  a heavily-doped  heterojunction 
transistor. 

5.5  Summary  and  Conclusions 

A theoretical  study  of  the  heavy  doping  effects  on  the  built-in 
potential,  depletion  layer  width,  current  density,  and  current  gain  for 
heavily-doped  abrupt  and  graded  hetero junction  devices  has  been  carried 
out  in  this  chapter.  A transition  region  is  modeled  to  extend  beyond 
the  space  charge  region  with  a constant  bandgap  gradient  for  the  graded 
het ero j u net i on  case.  Calculations  of  the  heavy  doping  effects  such  as 
carrier  degeneracy  and  bandgap  narrowing  were  carried  out  by  employing 
the  empirical  models  described  in  Section  2.2.  Numerical  results  for  a 
heavily-doped  A1  o.3Gao.7As_GaAs  heterojunction  transitor  are  given.  The 
results  of  the  present  study  show  that:  1)  bandgap  narrowing  and 
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carrier  degeneracy  have  important  effects  on  the  total  built-in 
potential  and  carrier  current  density,  as  well  as,  the  transistor 
current  gain,  8;  2)  the  effects  of  heavy  doping  on  the  depletion  layer 
width  and  depletion  capacitance  are  insignificant;  3)  the  effect  of 
carrier  degeneracy  is  the  dominant  heavy  doping  effect  in  the  heavily- 
doped  GaAs-Al GaAs  heterojunction  devices  due  to  the  very  small  electron 
effective  masses  in  GaAs  and  AlGaAs.  The  present  results  show  that  the 
new  calculated  values  for  the  built-in  potential,  carrier  current 
density,  and  transistor  current  gain  may  differ  significantly  from  those 
predicted  by  the  nondegenerate  theories. 
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p— GaAs  N-AJxGa^x  As 

DISTANCE  ► 

(NOT  DRAWN  TO  SCALE) 


Figure  5.1  Energy  band  diagram  before  formation  of  an  pn  GaAs- 
AlGaAs  heterojuction. 
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Figure  5.2  Schematic  representation  of  a GaAs-Al  Ga,  As  pngraded 
junction  in  thermal  equilibrium.  x "x 
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Figure  5.3 


Schematic  representation  of  a GaAs-Al  Ga, 
abrupt  heterojunction  in  thermal  equitibrT 
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Figure  5.4  Values  of  function  F(Y)  vs.  hole  or  electron  density  in 

the  heavily-doped  GaAs  at  T = 300  K.  Where  F(Y)  is  calcu- 
lated from  Eq. (2.3).  Solid  line  is  values  for  electrons, 
dash  line  for  holes. 
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Fig.  5 


.5  Values  of  function  F(Y)  vs.  hole  or  electron  density  in  the 
heavily-doped  A1  _ _Gan  ?As  at  T = 300  K.  Solid  line  is 
values  for  electrons , dash  line  for  holes. 
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5.6  Total  built-in  potential  vs.  emitter  dopant  density,  N , 
for  heavily-doped  abrupt  Al„  ~Ga„  ?As-GaAs  emitter-base 
hetero junction  at  T = 300  K.  Solid7 line  is  the  predicted 
values  by  using  Eq.  (5.19),  dash  line  is  calculated  by 
using  the  conventional  theory  from  Eq.  (5.18). 


Total  Built-in  Potential  ( V ) 


82 


Figure  5.7  Total  built-in  potential  vs.  base  dopant  density,  Nfl, 

for  heavily-doped  abrupt  Al_  -Gan  ,As-GaAs  emitter-base 
heterojunction  at  T = 300  K. 
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Emitter  Dopant  Density  ND(  cm'3) 


Figure  5.8  Depletion  layer  width  vs.  emitter  dopant  density,  N , 

for  heavily-doped  abrupt  Al„  ~Gan  ?As-GaAs  emitter-oase 
heterojunction  at  T = 300  K. 
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Figure  5.9  Hole  current  density  vs.  emitter  dopant  density,  Nn, 

for  heavily-doped  abrupt  A1  Gan  ?As-GaAs  emitter-base 
heterojunction  at  T = 300  K. 
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Figure  5.10  Current  gain  vs.  base  dopant  density,  N»,  for  an  npn 
transistor  with  wide  bandgap  AJ  ^Gag  yA§  emitter. 


CHAPTER  VI 

EFFECTS  OF  HEAVY  DOPING  ON  I-V  CHARACTERISTICS 
of  GaAs  JFET  AND  MESFET  DEVICES 


6.1  Introduction 

GaAs  JFETS  and  MESFETS,  particularly  the  E-JFET  and  E-MESFET 
devices,  show  great  promise  for  applications  in  high  speed,  low  power 
integrated  circuits  and  super  computers  [43,44].  In  addition,  for  some 
military  and  space  applications,  GaAs  ICs  are  particularly  important  for 
their  exceptional  radiation  hardness. 

To  increase  the  speed  of  FET  devices,  shortening  the  gate  length  is 
most  effective.  Many  experimental  results  [45-49]  show  a strong 
influence  of  short-channel  effects  if  the  channel  doping  is  less  than  1 
~ 2 x 10  cm”  and  the  gate  length  is  shortened  below  0.5  ym.  The  most 
pronounced  short  channel  effect  is  a strong  negative  shift  of  the 
threshold  voltage  Vy.  A shift  of  Vy  by  -1  to  -2  V are  reported  [48,49], 
respectively,  for  intermediate  doped  channels  (<  2 x 1017  cm"3)  when  the 
gate  length  is  decreased  from  more  than  1 wm  to  less  than  0.3  m while 
doping  concentrtion  and  channel  thickness  remain  constant. 

In  addition  to  the  lowering  of  the  threshold  voltage,  a slight 
decrease  of  transconductance,  and  the  lack  of  saturation  behavior  are 
also  observed  [50].  All  these  effects  have  to  be  taken  into  account  for 
device  layout  and  integrated  circuit  design.  One  effect,  for  instance, 
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is  that  devices  with  different  gate  lengths  need  variations  in  the 
production  processes  to  assure  equal  threshold  voltage.  On  the  other 
hand,  an  active  layer  thickness  laid  out  following  conventional  design 
rules  results  in  an  undesired  value  for  a submicrometer  gate's  threshold 
voltage  [51].  Therefore,  alleviation  of  these  effects  is  necessary. 

Theoretical  investigations  by  Drangeid  and  Sommerhalder  [52]  show 
that  the  performance  of  MESFET's  can  be  improved  by  reduction  of  the 
gate  length  if  the  dopiong  concentration  is  increased  correspondingly, 
so  that  the  ratio  of  gate  length  L to  the  channel  thickness  t is 
maintained  at  L/t  » 3.  Recently,  Daembkes  et  al . [51]  have  fabricated 
GaAs  MESFET  with  channel  dopant  density  up  to  5 x 1018  cm'3  along  with 
both  micrometer  and  submicrometer  gates.  They  reported  that  the  short- 
channel  devices  show  a good  saturation  behavior  and  no  shift  in  the 
threshold  voltage  compared  to  the  long-channel  devices.  Thus,  heavily- 
doped  FETs  demonstrate  a pronounced  alleviation  of  short-channel  effects 
as  experienced  for  1 x 10^  cm  3 doping  level.  Moreover,  it  was  shown 
that  the  use  of  heavily-doped  channels  enable  the  design  rules  for 
integrated  circuits  and  device  layout  to  be  maintained  down  to  less  than 
0.5  - ym  devices  [51]. 

Since  it  is  necessary  to  increase  the  dopant  density  in  order  to 
reduce  the  size  and  to  improve  the  electrical  characteristics  of  active 
devices  in  the  VLSI  (or  ULSI)  and  VHSIC  circuts,  it  became  more  and  more 
important  to  study  how  the  heavy  doping  effects  such  as  carrier 
degeneracy,  and  bandgap  narrowing  will  affect  the  performance  of  a 
heavily-doped  FET. 
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The  purpose  of  this  chapter  is'  to  theoretically  analyze  the  effect 
of  heavy  doping  on  the  threshold  voltage  and  I-V  characteristics  of  a 
heavily-doped  GaAS  JFET  and  MESFET.  Numerical  calculations  are 
performed  for  the  heavily-doped  E-JFET  and  E-MESFET  devices  using  a 
modified  expression  for  the  built-in  potential  described  in  Chapter 
II.  Empirical  approximation  for  the  reduced  Fermi -energy  is  presented 
by  [1],  and  the  GaAs  bandgap  narrowing  expression  is  given  by  [53],  The 
theoretical  treatment  which  includes  the  reformulation  of  the  built-in 
potential  for  a JFET  and  MESFET  are  given  in  Section  6.2.  The  result 
and  discussion  are  given  in  Section  6.3 

6.2  Theoretical  analysis 

The  disadvantage  of  the  short-channel  effects  of  a MESFET  can  be 
avoided  by  increasing  the  channel  doping  level  and  accordingly 
decreasing  the  channel  thickness  (i.e.,  the  ratio  of  gate  length  L to 
the  channel  thickness  t is  kept  at  L/t  >>  3).  Therefore,  the  current- 
voltage  relation  for  a uniformly-  and  heavily-doped  FET  may  still  be 
expressed  by  that  for  a long-channel  FET,  except  for  the  built-in 
potential  expression  which  should  be  replaced  by  Eq.  (2.12). 

Thus,  the  drain  current  for  a heavily-doped  short-channel  and  long- 
channel  FET  is  given  by 

■d  - 'p  (3VVp  - 2[<VD  - VG  + Vbi>3/2  - (vbl  - VG>3/2]/Vp/2l  (6-1) 

and  the  drain  saturation  current  for  a given  gate  voltage  is  obtained 
from  Eq.  (6.1)  as 
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DS  ’ V1  - 3»  + <VT  - V'V  + 2[1  + (VT  - VG)/Vp]3/2f  (6.2) 


2 2 3 

where,  Ip  (=  wpq  NQt  /6e$L)  is  the  pinch-off  current;  L is  the  channel 
length;  W is  the  channel  width;  y is  the  channel  mobility;  Vp  is  the 


channel  pinch-off  voltage;  NQ  is  the  channel  doping  density  for  the 
N-channel  FET;  VQ  is  the  applied  drain  voltage;  e$  is  the  material 
dielectric  constant;  t is  the  channel  thickness,  and 


is  the  threshold  voltage.  The  built-in  potential,  Vbi , of  the  gate  and 
channel  junction  in  a conventional  JFET  device  is  given  by  [54] 


where  Nc  is  the  effective  density  of  conduction  band  states,  Egi  is  the 
intrinsic  bandgap  of  semiconductor,  and  Vth  (=  kT/q)  is  the  thermal 
voltage.  Since  the  junction  between  gate  and  channel  in  a JFET  is  a 
heavily-doped  pn  junction,  the  built-in  potential,  VbiJ,  for  a heavily- 
doped  n-channel  JFET  should  be  modified  to  include  the  heavy-doping 
effects  as  described  in  Chapter  II.  It  can  be  rewritten  here  as: 


(6.4) 


(6.5) 


and 


(6.6) 
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where  Yp  = ND/NC,  Yp  = P/Nv,  P is  the  hole  density  in  the  heavily-doped 

p-gate  region,  and  F(Yn)  (or  (F(Yp))  can  be  approximated  by  Eq. 

(2.13).  Fig.  5.4  shows  the  calculated  values  of  F(Yn)  and  F(Y  ) for 

GaAs,  values  of  F(Yn)  vary  from  0.075  to  6.14  for  electron  densities  in 

the  range  of  1 x 1017  to  1 x 1019  cm-3,  and  values  of  F(Y  ) vary  from 

0.05  to  4.34  for  hole  densities  in  the  range  of  1 x 1018  cm"3  to  1 x 
20  -3 

10  cm"  . AEgN  and  AEgp  are  the  bandgap  narrowing  in  the  n-  and  p-type 
semiconductor,  respectively  [53].  Since  the  effective  mass  of  electron 
in  GaAs  is  very  small,  the  effect  of  bandgap  narrowing  is  negligible  for 
the  present  case.  On  the  other  hand,  carrier  degeneracy  is  the  dominant 
heavy  doping  effect  in  a heavily-doped  n-channel  GaAs  FET. 

The  built-in  potential,  VbiM,  of  the  Schottky  barrier  in  a heavily- 
doped  n-channel  MESFET  can  be  expressed  by  [54] 

VbiM  = VBN  + EFf/q  (6*6) 

where  the  bandgap  narrowing  is  very  small  and  can  be  neglected.  VBN  is 
the  barrier  height  of  the  metal -semiconductor  contact.  The  Fermi -energy 
for  an  n-type  nondegenerate  semiconductor  is  given  by  EpN  = kT  An(n/Nc), 
and  for  a degenerate  semiconductor,  it  can  be  approximated  by  [1] 

ND  Nd  . 

EFN  ~ kT[£n(jj— ) + F(-^— )]  (6.7) 

c c 

Thus,  the  transconductance  of  an  FET  in  the  saturation  region  can  also 
be  expressed  by 
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V ' W1  - [<VbiJ(M)  - Vg’/Vp11/2}  <6-8) 

where  VbiJ  and  VbiM  are  the  built-in  potential  given  by  Eqs.  (6.5)  and 
(6.6),  respectively;  and 

gmax  = gNDtyw/L  • (6.9) 

p is  dependent  on  channel  doping  density  and  can  be  described  by 
Hilsum's  formula  [55]. 

6.3  Results  and  Discussion 

Calculations  of  the  effects  of  heavy  doping  on  the  threshold 
voltage  shift,  I-V  characteristics,  and  transconductance  of  a JFET  and 
MESFET  are  carried  out  at  T = 300  K by  using  equations  described  in 
Section  6.2.  Figs.  6.1  and  6.2  show  the  effects  of  carrier  degeneracy 
and  bandgap  narrowing  on  the  threshold  voltage  of  a heavily-doped 
n-channel  GaAs  E-JFETS  and  E-MESFETS.  Fig.  6.1  shows  the  threshold 
voltage  shift  vs.  p-gate  region  density,  NA,  for  an  n-channel  E-JFETS, 
in  which  channel  dopant  density,  NQ,  is  assumed  to  be  2 x 1018  cm-3,  and 
the  threshold  voltage  is  about  0.5  volt.  The  threshold  voltage  shift 
vs.  channel  dopant  density  for  n-channel  GaAs  E-MESFETS  is  shown  in  Fig. 
6.2,  where  the  barrier  height,  VBN,  is  assumed  to  be  0.8  volt,  and  the 
threshold  voltage  is  about  0.3  volt  for  E-MESFETS.  From  Figs.  6.1  and 
6.2,  it  is  clearly  show  that  the  effects  of  heavy  doping  tend  to 
increase  the  threshold  voltage  as  high  as  36%  and  80%  for  E-JFETS  and 
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E-MESFETS , respectively.  On  the  contrary,  the  short-channel  effects 
tend  to  reduce  the  threshold  voltages  for  short-channel  FETS.  Thus, 
raising  the  channel  doping  level  for  short-channel  FETs  not  only  can 
avoid  the  disadvantages  of  the  short-channel  effects  of  a FET  and 
satisfy  the  criteria  of  L/a  » 3,  but  also  can  act  as  an  adjustment  of 
the  FET  threshold  voltage.  The  I-V  characteristics  of  an  n-channel 
E-MESFET  with  gate  voltage  as  a parameter  is  shown  in  Fig.  6.3.  Where 
channel  dopant  density,  ND,  is  assumed  to  be  1.5  x 10"18  cm"8,  channel 
depth  a = 225  A,  and  barrier  height,  VBN,  is  assumed  to  be  0.8  volt. 
Solid  lines  are  the  predicted  values  by  Eqs.  (6.1),  (6.2),  and  (6.7). 
Dash  lines  are  the  values  calculated  by  the  conventional  nondegenerate 
theory  from  [54],  The  results  show  that  the  drain  current  of  E-MESFETS 
will  decrease  significantly  when  carrier  degeneracy  effect  is  included 
in  the  calculations. 

Table  1 shows  the  calculated  results  of  drain  saturation  current 
and  saturation  transconductance  for  several  GaAs  D-MESFETS  and  E-MESFETS 
with  specific  device  parameters  and  gate  voltages.  It  is  noted  that  the 
effect  of  carrier  degeneracy  plays  an  important  role  on  the  drain 
saturation  current,  Iq<j,  and  saturation  transconductance,  gms,  of 
D-MESFETS  as  well  as  E-MESFETS  via  the  built-in  potential  presented  in 
the  heavily-doped  junction  between  gate  and  channel.  It  is  expected 
that  the  carrier  degeneracy  will  also  have  significant  influence  on  the 
IpB  and  gms  of  E-JFETS.  On  the  contrary,  carrier  degeneracy  will  have 
only  little  influence  on  the  D-JFETS. 
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6.4  Summary  and  discussion 

In  short,  the  considerations  of  short-channel  effects  on  the 
heavily-doped  FETs  are  justified  by  Daembkes  et  al . [54]  and  were  shown 
to  be  not  important  for  heavily-doped  short-channel  FETS.  The  effects 
of  heavy  doping  on  the  I-V  characteristics  of  JFETS  and  MESFETS  have 
been  theoretically  investigated  in  this  chapter.  It  is  shown  that 
carrier  degeneracy  has  a significant  effect  on  the  threshold  voltage  and 
the  I-V  characteristics  of  a heavily-doped  GaAs  E-JFETS  and  E-MESFETS, 
while  the  bandgap  narrowing  effect  has  negligible  influence  on  the 
electrical  characteristics  because  of  the  very  small  electron  effective 
mass. 


94 


TABLE  6.1 


Comparisons  of  calculated  drain  saturation  current  and  saturation 
transconductance  of  GaAs  JFETs  and  MESFETs  with  specified  device 
parameters  and  applied  gate  voltage. 


Device  E-MESFET 

E-MESFET 

D-MESFET 

D-MESFET 

n-channel 

n-channel 

n-channel 

n-channel 

ND=1.5xl018cirT3 

ND=5xl018cm-3 

Nq=1 .5xl0-*-8cm~3 

ND=5xl018cm"3 

#Vbn  = 0.8 

VBN  = 0,8 

VBN  = 0,8 

VBN  " °*8 

a = 225  A 

a = 125  A 

a = 405  A 

a = 225  A 

Vq  = 0.7  volt 

Vq  = 0.7  volt 

Vq  = 0.0  volt 

Vq  = 0.0  volt 

(Theory) 

parameter 

(Degenerate) 

IDS/Ip  0.47 

0.285 

0.211 

0.182 

(Nondegenerate) 

W!P  °-7 

0.823 

0.251 

0.281 

(Degenerate) 

0ms70max  0.481 

0.353 

0.297 

0.272 

(Nondegenerate) 

Oms^max  0.634 

0.731 

0.327 

0.349 

# in  unit  of  volt 


Threshold  Voltage  Shift  aVt  ( Volt  ) 
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Figure  6.1  Threshold  voltage  shift  vs.  p-gate  dopant  density 

for  n-channel  E-JFETs,  in  which,  channel  dopant  density, 
N , is  assumed  to  be  2 x 10^8  cm~3. 


Threshold  Voltage  Shift  AVt  ( Volt  ) 
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Figure  6.2  Threshold  voltage  shift  vs.  channel  dopant  density  for  n- 
channel  E-MESFETs. 


Drain  Current  (ln/l_)  (in  unit  of  I 
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Figure  6.3  I-V  characteristics  of  an  n-channel  E-MESFET  with  gate 
voltage  as  a parameter.  Where,  channel  dopant  density 
Ng  = 1.5  x ICr^  cm"3,  chnnel  depth  a = 225  A,  and  ba- 
rrier height  Vg^  is  assumed  to  be  0.8  volt. 


CHAPTER  VII 

EXPERIMENTAL  METHODS  FOR  STUDYING  THE 
METALLIC  CONTAMINATION  IN  SILICON  DEVICES 


Experimental  tools  employed  in  the  study  of  the  effects  of  metallic 
contamination  in  electrical  characteristics  of  semiconductor  devices 
include  the  current-voltage  (I-V)  measurement,  capacitance-voltage  (C-V) 
measurement,  thermally  stimulated  capacitance  (TSCAP)  measurement,  and 
deep  level  transient  spectroscopy  (DLTS)  measurement.  From  these 
measurements,  one  can  determine  the  diode  characteristics  and  the  defect 
parameters  such  as  defect  energy  levels,  defect  densities,  and  capture 
cross  sections.  Each  of  these  measurement  techniques  is  described  in 
the  following  sections. 

7.1  Current-voltage  (I-V)  measurement 

Measurement  of  the  current-voltage  (I-V)  characteristics  under 
forward  bias  condition  can  yield  useful  information  concerning  the 
conduction  mechanisms,  recombination  processes  in  the  space  charge 
region  of  a p-n  junction  diode  or  a Schottky  barrier  diode.  For  a good 
p-n  diode  with  no  surface  leakage,  the  total  current  is  composed  of  the 
diffusion  current  in  the  quasi -neutral  region  (QNR)  and  the 
recombination  current  in  the  space  charge  (SCR).  When  bulk  diffusion 
current  component  dominates,  the  current  expression  is  given  by  [9] 
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If  = Id[exp(Va/VT)  - 1]  (7.1) 

Where  Va  is  the  applied  voltage,  VT  = kT/q,  and  Id  is  the  magnitude  of 
the  saturation  diffusion  current. 

'd  = <"’iA[<Dn/LnNA>  * <VLPND)]  <7-2> 

where  ni  is  the  intrinsic  carrier  density;  A is  the  diode  area;  Dn  (Dp) 
is  the  diffusion  constant  for  electrons  (holes);  Lp  (Lp)  is  the 
diffusion  length  for  electrons  (holes);  and  NA  (ND)  is  the  acceptor 
(donor)  density.  If  bulk  generation-recombination  current  dominates, 
then  the  current  expressin  is  given  by  [9]: 

lr  = Trg  exP  (y2vT)  (7.3) 

where  Ipg  is  the  magnitude  of  the  generation-recombination  current. 


I 


rg 


qn.W/2xe 


(7.4) 


In  Eq.  (7.4),  W denotes  the  depletion  layer  width;  = 
effective  carrier  lifetime  in  the  space  charge  region 
the  lifetime  of  electrons  (holes),  defined  by: 


(TnTp)1/2  is  the 

[56];  xp  (xp)  is 


x 

n 


1/(NtVvth>) 


(7.5) 
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Where  Nt  is  the  trap  density,  (d  ) is  the  capture  cross  section  for 
electrons  (holes);  <vt^>  is  tHe  average  thermal  velocity. 

The  total  current  can  be  expressed  as: 

It  = If  + Ip  » IQ  exp  (Va/n*VT)  (7.6) 

Where  IQ  is  the  saturation  current;  n*  is  the  ideality  factor.  The 
diode  ideality  factor,  n can  be  used  to  identify  the  dominant  current 
component  in  a p-n  diode.  Inspection  of  Eqs.  (7.1)  and  (7.3)  show  that 
the  bulk  diffusion  current  depends  more  strongly  on  temperature  than  the 
recombination  current  in  the  junction  SCR.  Since  the  recombination 
current  in  the  junction  SCR  is  inversely  proportional  to  the  effective 
carrier  lifetimes  (and  hence  directly  related  to  the  defect  density  in 
the  transition  region),  measurements  of  I-V  characteristics  under 
forward  bias  conditions  would  allow  us  to  determine  the  effective 
1 ifetimes  in  Si . 

7 .2  Capacitance-voltage  (C-V)  measurement: 

The  capacitance-voltage  (C-V)  measurement  can  be  used  to  determine 
the  background  doping  concentration  in  the  n-  or  p-  Si  epitaxial  layers 
using  a Schottky  barrier  structure  or  a one-sided  abrupt  p+-n  (or  n+-p) 
junction.  The  deplation  capacitance  across  the  Schottky  barrier  diode 
is  given  by: 


1/2 


C(VT)  = e$A/W  = A{qesND/[2(4.i  + Vp  - kT/q)  ]} 


(7.7) 
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where  eg  is  the  dielectric  constant  of  semiconductor,  ^ is  the  built-in 
potential,  and  Vp  is  the  applied  reverse  voltage.  Eq.  (7.7)  shows  that 
the  depletion  capacitance  of  a Schottky  diode  is  proportional  to  the 
square  root  of  dopant  concentration  and  inversely  proportional  to  the 
square  root  of  the  applied  voltage.  If  the  inverse  of  the  capacitance 

O 

square  (C_<:)  is  plotted  as  a function  of  the  reverse  bias  voltage  Vp, 

then  the  background  concentration  can  be  calculated  from  the  slope  of 

-2 

C vs.  Vp  using  the  following  expression: 

C“2(Vr)  = [2/ (qG$ A2Nd ) ](*.  + Vp)  (7.8) 

The  intercept  of  C ^ vs.  Vp  plot  with  the  voltage  axis  yields  values  of 
4> -j  which  is  related  to  the  barrier  height  of  a Schottky  diode  by: 

*BN  = +1  + Vn  + IT  - A<f>  (7-9) 

where 

Vn  = Ec  " (kT/q)  in  (Nd/Nc)  (7.10) 

and  A4>  is  the  image  lowering  potential  of  a Schottky  diode;  Nc  is  the 
density  of  state  in  the  conduction  band. 

7.3  Thermally  stimulated  capacitance  (TSCAP)  measurement: 

Thermally  stimulated  capacitance  measurements  (TSCAP)  add  another 
improvement  over  the  stepped  C-V  measurement  technique  as  reported  by 
Sah,  Senechal , and  Basinski  [57,58].  The  main  theoretical  aspect  which 
controls  the  TSCAP  results  [59-61]  is  given  by 
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el  = f (E,T  ,t)  = 


-(Ec  - Et)/kT 


(7.11) 


where  ejj  is  the  electron  thermal  emission  rate,  Ec 


E^.  is  the  trap 


depth  measured  from  the  band  edge,  g is  a degeneracy  factor  of  the 
defect  level,  k is  the  Boltzmann  constant,  <vth>  is  the  carrier  thermal 
velocity,  Nc  is  the  effective  density  of  conduction  band  states,  T is 
the  absolute  temperature,  and  an  is  the  thermal  capture  cross-section 
for  electrons.  If  we  ignore  the  temperature  dependence  of  some  of  the 
terms  given  in  eq.  (7.11),  the  general  expression  indicates  that  each 
trap  will  have  a unique  value  of  (Ec  - E^.)  which  will  make  the  exponen- 
tial term  of  eq.  (7.11)  maximum  at  a corresponding  unique  temperature. 
This  means  that  when  the  unique  temperature  for  each  trap  is  reached, 
the  emission  for  that  trap  will  increase  exponentially.  However,  the 
source  of  the  emission  is  limited  and  must  follow  a certain  time  cons- 
tant, and  a peaked  emission  profile  of  the  trap  will  result  (as  shown  in 
Fig.  7.1)  for  a linearly  increasing  temperature.  Fig.  7.2  shows  a 
typical  TSCAP  scan  that  shows  only  majority  carrier  traps  for  a Schottky 
barrier  (three  single  trap  changes  and  a combined  unknown  at  77°K).  A 
derivative  (dC/dT)  of  Fig.  7.2  will  then  result  in  a capacitance  versus 
temperature  plot  similar  to  Figure  7.2.  It  is  also  possible  that  some 
shallow  traps  may  be  present  that  would  emit  electrons  at  temperatures 
lowr  than  77K.  In  this  case,  the  net  trap  density  can  be  determined  by 
recording  the  change  (AC  = 0.7  pf  in  Fig.  7.2)  in  the  capacitance 
between  t = 0+  (immediately  after  reverse-bias)  and  at  t = °°  (after  it 
ceases  to  change),  while  holding  the  temperature  at  77  K. 
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The  amplitude  of  the  capacitance  step  (as  shown  in  Fig.  7.2)  is 
directly  proportional  to  the  trap  density.  The  trap  density  for  n-Si 
can  be  calculated  from  the  following  expression: 

Nt  = V2AC/Co}  (7.12) 

where  CQ  is  the  depletion  layer  capacitance,  and  AC  is  the  capacitance 
change  due  to  the  majority  or  minority  carrier  emission.  Thus,  knowing 
Nq  (or  and  CQ  at  the  temperature  where  the  capacitance  step  was 
observed,  the  trap  density  can  be  calculated  from  Eq.  (7.12).  Note  that 
Eq.  (7.12)  is  valid  only  for  the  case  when  Nt  < 0.1  NQ.  For  the  case  of 
large  trap  density  with  > 0.1  Nq,  a more  exact  expression  should  be 
used  instead. 

7 .4  Deep  level  transient  spectroscopy  (DLTS ) : 

The  Deep-Level -Transient-Spectroscopy  (DLTS)  experiment  is  a high 
frequency  (20  MHZ)  transient  capacitance  technique,  which  was  introduced 
first  by  Lang  in  1974  [62,63].  The  DLTS  scan  displays  the  spectrum  of 
deep  level  traps  in  the  forbidden  gap  of  a semiconductor  as  positive  or 
negative  peaks  on  a flat  baseline  as  a function  of  temperature. 
Although  this  kind  of  measurement  is  time  consuming,  it  offers  several 
advantages.  It  is  sensitive,  easy  to  analyze  and  capable  of  measuring 
the  traps  over  a wide  range  of  depth  in  the  forbidden  gap. 

Fig.  7.4.  shows  the  schematic  block  diagram  of  a DLTS  system  used 
in  this  study.  The  DLTS  system  consists  of  a sensitive  capacitance 
measurement  apparatus  with  good  transient  response,  a pulse  generator,  a 
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boxcar  with  dual -gated  signal  integrator,  an  X-Y  recorder,  and 
temperature  control  system  for  heating  and  cooling.  By  properly 
changing  the  experimental  parameters  it  is  possible  to  measure  the 
following  parameters: 

- The  minority  and  majority  carrier  traps. 

- The  activation  energy  of  each  defect  level. 

- The  defect  concentration  which  is  directly  proportional 
to  the  peak  height. 

- The  defect  concentration  profile. 

- The  electron  and  hole  capture  cross-sections  for  each 
trap. 

7.4.1  Principles  of  the  DLTS  technique 

The  capacitance  transient  is  associated  with  the  return  to  thermal 
equilibrium  of  the  carrier  occupancy  in  a trap  level  following  an 
initial  nonequilibrium  condition.  The  polarity  of  the  DLTS  peak  depends 
on  the  capacitance  change  after  trapping  the  minority  or  majority 
carriers.  Because  an  increase  in  trapped  minority  carriers  in  the  SCR 
would  result  in  an  increase  in  the  junction  capacitance,  the  minority 
carriers  trapping  will  produce  a positive  polarity  peak,  and  vice 
versa.  For  example,  in  a p+n  junction  diode,  the  SCR  extends  mainly 
into  the  n-region,  and  the  local  charges  are  due  to  positively  charged 
ionized  donors.  If  a forward  bias  is  applied,  the  minority  carriers 
(holes)  will  be  injected  into  this  region.  Once  the  holes  are  trapped 
in  a defect  level,  the  net  positive  charges  in  such  region  will 
increase.  This  results  in  a narrower  SCR  width  which  implies  a positive 
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capacitance  change.  Thus,  the  DLTS  signal  will  have  a positive  peak. 

Similarly,  if  the  majority  carriers  are  injected  into  this  region  and 
captured  by  the  majority  carrier  traps,  which  reduces  the  local  charges, 
the  SCR  width  will  be  wider,  implying  a decrease  of  the  junction 
capacitance.  Therefore,  the  majority  carriers  trapping  will  result  in  a 
nagative  DLTS  peak.  The  same  argument  can  be  applied  to  the  n+p 
junction  diodes.  All  of  the  samples  used  in  the  DLTS  measurements  are 
p+n  diodes  so  that  a positive  peak  represents  the  hole  trap,  and  a 
negative  peak  represents  the  electron  trap. 

7.4.2  Minority  carrier  injection 

Fig.  7.5  shows  the  injection  of  minority  carriers.  Fig.  7.5(a)  is 
the  injection  pulse  where  V = Vp  > 0 during  the  pulse  (ta  < t < tb)  and 
V = VR  < 0 outside  the  time  interval  (ta,tb).  Fig.  7.5(b)  is  a 
simplified  energy  band  diagram  (the  band  bending  due  to  the  junction 

electric  field  is  omitted)  in  quiescent  reverse  bias  condition  (t  < 

ta)»  Fig.  7.5(c)  is  the  saturation  injection  pulse  (ta  < t < tb). 
During  the  pulse,  the  minority  carrier  traps  Et  are  filled  by  holes. 

Fig.  7.5(d)  is  the  transient  process  during  the  period  (t  > tb);  the 
trapped  hoes  begin  to  emit  from  the  trap  centers  to  valence  band  and 
then  swept  out  of  the  SCR  by  the  built-in  electric  field  that  causes  the 
capacitance  transient. 

7.4.3  Majority  carrier  injection 

Fig.  7.6  shows  the  injection  of  majority  carriers.  Fig.  7.6(a)  is 
the  injection  pulse  where  V = Vp  = 0 during  the  pulse  (tg  < t < tb)  and 
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V = VR  < 0 outside  the  time  interval  (ta,  tb).  Fig.  7.6(b)  is  a 
simplified  band  diagram  in  quiescent  reverse  bias  condition  (t  < ta). 
Fig.  7.6(c)  is  the  majority  carrier  injection.  During  the  pulse  period 
^a  t < tb^ » the  majority  carriers  were  injected  into  the  SCR  and 
captured  by  the  majority  trap.  Fig.  7.6(d)  is  the  transient  process  for 
the  period  (t  > tb)  in  which  the  captured  electrons  begin  to  emit  from 
the  trap  centers  to  conduction  band  and  then  swept  out  of  the  SCR  by  the 
built-in  electric  field,  resulting  in  the  observed  capacitance 
transient. 

7.4.4  Trap  concentration 

The  defect  concentration  is  directly  proportional  to  the  peak 
height  as  described  before,  and  the  peak  height  is  proportional  to  the 
capacitance  change  AC(0)  [AC(0)  is  shown  in  Figs.  7.5(d)  and  7.6(d)]. 
Therefore,  the  defect  concentration  Nt  is  proportional  to  AC(0)  to  be 
derived  as  follows:  Let  C(t)  be  the  capacitance  transient  [64,65]  which 

is  proportional  to  the  electron  (hole)  emitted  to  the  conduction 
(valence)  band.  Then 

C(t)  = A{qes[ND  - Nt  exp  (-t/x) ]/ [2( + Vp  + kT/q)  ] }1/2 

= CQ{1  - [Nt  exp  ( -t/x)/ ND ] }1/2  (7.13) 

where  t is  time;  x is  the  carrier  emission  time  constant;  CQ  = C(V  ),  as 
shown  in  Eq.  (7.7),  is  the  junction  capacitance  at  the  quiescent  reverse 
bias  condition.  For  N^/N^  <<  1,  Eq.  (7.13)  reduces  to  a simple  form  as: 
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C(t)  = C0[l  - Nt  exp  (-t/T)/2ND] 


(7.14) 


Eq.  (7.14)  can  be  rewritten  as 


Nt  exp  (-t/x)  = [2AC(t )/CQ ] • Nd 


(7.15) 


where  AC(t)  = CQ  - C(t).  From  the  DLTS  technique,  AC(0)  can  be 
measured.  The  junction  capacitance  CQ  and  the  background  concentration 
can  be  obtained  from  C-V  measurements.  Thus,  the  defect  concentration 
Nt  can  be  calculated  easily  by  using  Eq.  (7.15)  at  t = 0. 

7.4.5  Activation  energy  of  a deep  level  trap 

As  shown  in  Fig.  7.7,  AC(t)  begins  to  decay  after  the  injection 
pulse  is  over.  The  decay  associates  with  a specific  time  constant  which 
is  equal  to  the  reciprocal  of  the  emission  rate.  For  an  electron  trap, 
the  emission  rate  "en"  is  a function  of  temperature,  coefficient  and 
activation  energy,  and  can  be  expressed  by  [11,66] 


where  Et  is  the  activation  energy  of  the  trap;  g is  the  degeneracy 
factor,  0^  is  the  electron  capture  cross-section  which  is  temperature 
dependent  [67] 


e 


n 


= (Vvth>Nc/g)  * exp  t(Ec  " Et)/kTJ 


(7.16) 


°n  = a«  ’ e*P  (-AEb/kT) 


(7.17) 
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where  a„  is  the  capture  cross-section  at  very  high  temperature;  AEb  is 
activation  energy  of  the  capture  cross-section  for  the  trap.  en  can  be 
written  as 


en  = BT2  • exp  {[Ec  - (Et  + AEb) ]/kT} 
= BT2  • exp  [(Ec  - Em)/kT] 


(7.18) 


where  B is  the  proportionality  constant  and  is  independent  of 
temperature.  From  this  relation,  en  increases  with  increasing 
temperature.  The  capacitance  transient  is  rearranged  from  Eq.  (7.14)  as 


C(t)  = Co(Nt/2N0)  • exp  (-t/x) 
= AC(0)  • exp  (-t/x) 


(7.19) 


where  x = en_1.  Eq.  (7.19)  implies  a faster  decay  of  AC(t),  as  is  shown 
in  Fig.  7.6,  at  high  temperature  and  a smaller  decay  time  constant  which 
implies  a faster  decay  of  AC(t)  as  shown  in  Fig.  7.6.  The  same  argument 
can  also  be  applied  to  the  hole  trap.  The  only  difference  is  that  AC(t) 
has  an  opposite  polarity. 

The  following  is  the  procedures  to  deduce  the  activation  energy  of 
a defect  level.  We  first  set  t^  and  t2  in  a dual -gated  integrator 
boxcar,  then  we  write 
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0(1^)  = C(0)  • exp  (-tj/x)'  (7.20) 

C(t2)  = C(0)  • exp  (-t2/x)  (7.21) 

The  DLTS  scan  along  the  temperature  axis  is  obtained  by  taking  the 
difference  of  Eq.  (7.20)  and  (7.21)  which  yeilds: 

S(T)  = AC(0)[exp(-t1/x)  - exp(-t2/x)]  (7.22) 

The  maximum  emission  rate,  x"*x,  is  obtained  by  differentiating  S(T) 
with  respect  to  x,  and  let  it  equal  zero  [dS(T)/dx  = 0].  The  result 
yields: 


Tmax  = <‘l  - t2)/*n(t1/t2)  (7.23) 

Under  this  condition,  S(T)  reaches  its  maximum  value  at  a specific 
temperature.  The  emission  rate  is  given  by  ep  = l/xmax  for  each  t1  and 
t2  setting.  By  changing  tj  and  t2  several  times,  a set  of  temperatures 
corresponds  to  this  set  of  xmgx  (or  emission  rate  en)  can  be  obtained. 
From  the  Arrhenius  plot  [i.e.,  £n(en/T2)  vs.  1/kT],  the  activation 
energy  of  the  trap  can  be  calculated  from  the  slope  of  this  plot. 
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Figure  7.1  Rate  of  ionized  traps  vs.  a linear  temperature 
increase. 

C(T) 


Figure  7.2  Typical  TSCAP  linear  temperature  scan  of  a Schottky 
diode. 
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Figure  7.3  dC/dT  (derivative)  vs.  temperature  of  Figure  7.2 
showing  the  three  trap  peaks  to  give  additional 
visibility  and  resolution  over  the  TSCAP  plot. 
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Figure  7.4  The  block  diagram  of  the  DLTS. 
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Figure  7.5  The  process  of  minority  carrier  injection. 
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Figure  7.6  The  process  of  majority  carrier  injection. 
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Figure  7. 


AC<t) 


7 The  transient  capacitance  decay  exponentially  after  the 
injection  pulse  is  over.  The  higher  temperature  corres- 
ponds to  the  smaller  decay  time  constant. 


CHAPTER  VIII 

THE  EFFECTS  OF  METALLIC  CONTAMINATION 
ON  THE  DEVICE  PERFORMANCE 


8.1  Introduction 

In  the  IC  fabrication,  the  unintentional  metallic  contamination  may 
introduce  into  silicon  substrate  during  high  temperature  processes. 
Most  of  these  metallic  contamination  is  associated  with  the  transition 
metals. 

The  3d  transition  metals  such  as  Ti  and  Ni , and  the  noble  metals 
such  as  Cu,  Ag,  and  Au  will  diffuse  interstiti al ly  into  silicon  and  stay 
predominantly  in  these  sites  during  high  temperature  processes.  These 
impurities  may  introduce  deep  donor  level  or  deep  acceptor  centers  in 
silicon  [68].  In  bipolar  devices,  these  deep  levels  may  act  as 
recombination  centers  under  forward  bias  condition.  The  effect  would  be 
to  reduce  the  current  gain  of  a transistor,  or  to  decrease  its  switching 
time.  In  a reverse  biased  p-n  junction,  these  deep  level  traps  act  as 
current  generation  centers,  and  increase  reverse  leakage  and  noise  in  a 
p-n  diode,  or  a transistor.  Thus,  the  knowledge  about  the  behavior  and 
the  influence  of  various  transition  metal  impurities  in  a semiconductor, 
especially  in  silicon,  is  very  important  for  improving  the  device 
performance. 

The  process-related  defects  in  silicon  has  been  actively 
investigated  by  many  researchers  [68-72].  Benton  and  Kimerling  [71]  use 
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the  capacitance  transient  spectroscopy  monitors  electrically  active 
defects  in  the  depletion  region  of  a p-n  junction.  They  recorded  the 
capacitance  transient  spectra  for  metallic  impurities  in  silicon  and 
correlated  it  with  sample  spectra  to  identify  trace  contamination  in  the 
tested  devices.  The  DLTS  exhibits  a characteristic  thermal  scanning 
peak  for  each  defect.  A complete  DLTS  analysis  determines  the  defect 
state  concentration,  the  activation  energy  for  the  trap  levels,  the 
carrier  capture  cross  section,  and  the  depth  distribution  of  the  defect 
in  a semiconductor.  Values  of  of  activation  energy  Ey  and  majority 

carrier  capture  cross  section  (a)  give  a unique  set  of  parameters  for 
each  defect.  Thus,  from  the  DLTS  spectra  and  the  measured  defect 
parameters  one  can  determine  the  specific  impurities  present  in  each 
device  [71]. 

In  this  chapter,  studies  of  metallic  contamination  in  low  noise 
silicon  JFET's  are  presented.  Experimental  tools  employed  in  this  study 
include  the  current  voltage  (I-V),  capacitance  voltage  (C-V),  thermally 
stimulated  capacitance  (TSCAP)  and  deep-level  transient  spectroscopy 
(DLTS)  measurements.  The  correlations  among  the  trap  density,  the 
reverse  leakage  current,  and  the  1/f  noise  are  also  investigated. 

The  device  structures  used  in  this  study  are  described  in  section 
8.2.  Experimental  results  obtained  from  the  I-V,  C-V,  TSCAP,  and  DLTS 
measurements  on  vertical  P/N  and  N/P  diodes  as  well  as  the  MOS 
capacitors  built  on  N-  and  P-collectors  are  presented  and  discussed  in 
section  8.3.  Section  8.4  gives  the  summary  and  conclusions. 
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8.2  Tested  devices  and  structures 

The  large  area  junction  diodes  and  MOS  capacitors  used  in  ths  study 
are  fabricated  from  ion-implanted  JFETs  process  (run  #JFET2-1358)  and 
found  on  the  test  mask  set  shown  in  Figure  8.1.  The  large  area  junction 
diodes  are  illustrated  in  Figure  8.2  and  8.3.  Figure  8.2  is  a large  PN 
junction  structure,  either  N diffused  into  P or  P diffused  into  N,  with 
a maximum  bottom  plate  area  and  a minimum  sidewall  area.  The  structure 
shown  in  Figure  8.3  differs  from  Figure  8.2  in  that  it  has  a maximum 
sidewall  area  and  a minimum  bottom  plate  area.  Figure  8.2  is  referred 
to  as  a bottom  plate  (BP)  diode  structure,  while  Figure  8.3  is  referred 
to  as  a sidewall  (SW)  diode  structure. 

Table  8.1  summarized  the  tested  devices  showing  their  structures 
and  applications.  Figure  8.4  illustrates  the  cross-sectional  view  of 
the  four  diode  structures  Dl,  D2,  D3,  and  D4.  Two  MOS  capacitor 
structures  are  shown  in  Figure  8.5. 

Several  die  from  wafers  of  the  JFET2-1358  run  were  characterized  in 
this  study.  The  group  B diodes  from  die  B-5  and  B-8,  are  bonded  out 
from  the  N+/Pb  BP,  N+/Pb  SW  and  P+/Nb  BP  diodes.  The  group  C diodes 
from  die  C-9  and  C-10  are  bonded  out  from  the  Nb/Pc  BP,  Nb/Pc  SW  and 
Pb/N-  and  Pc  bottom  plates.  The  P2P1  parts  are  bonded  out  from  the 
Nb/Pc  BP  and  Pb/N-  BP  diodes. 

8.3  Experimental  observation 

In  this  section,  the  experimental  results  obtained  from  the 
metallic  contamination  study  in  the  devices  from  JFETS  runs  are 
discussed. 
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8.3.1  Results  and  discussions  on  group  B (P1PB)  diodes 

The  forward  I-V  characteristics  of  both  B-5  and  B-8  N+/ Pb  BP 
diodes,  N+/Pb  SW  diodes,  P+/Nb  SW  and  P+/Nb  BP  diodes  are  shown  in 
Figures  8.6  and  8.7  with  good  diode  ideality  factors  ranging  from  1.03 
to  1.12.  However,  the  reverse  I-V  characteristic  curves,  see  Figure 
8.8,  show  a large  variation  in  leakage  current,  as  much  as  five  orders 
of  magnitude,  between  diodes  made  from  the  same  diffusion  on  two 
different  die.  Figure  8.8  shows  that  the  N+/Pb  SW  and  the  N+/Pb  BP 
diodes  from  the  B-5  group  have  the  highest  leakage,  followed  by  the 
P+/Nb  BP  diodes  from  the  group  B-5  and  B-8  diodes.  The  N+/Pb  SW  and  BP 
diodes  from  the  B-8  die  have  the  lowest  reverse  leakage  current. 

C-V  measurements  were  made  on  N+/Pb  BP  and  SW  diodes  and  P+/Nb  BP 
and  SW  diodes,  and  the  results  are  shown  in  Figures  8.9  and  8.10.  The 
dopant  density  in  the  base  region  of  these  diodes  was  determined  from 
1/C  vs.  VR  plots  and  the  results  are  summarized  in  Table  8.2  and  8.3. 
The  dopant  density  was  found  to  vary  from  2 x 10^  cm"3  in  the  Nc  region 
to  5 x 1017  cm"3  in  the  Nb  region. 

Deep-Level  Transient  Spectroscopy  (DLTS)  measurements  on  several 
group  of  diodes  were  carried  out.  The  N+/Pb  BP  and  SW  diodes  on  the  B-8 
die  had  relatively  low  reverse  leakage  current  as  compared  to  the  same 
diodes  on  the  B-5  die  as  is  seen  in  Figure  8.8.  The  higher  reverse 
leakage  current  observed  in  the  B-5  N+/Pb  BP  and  SW  diodes  was 
attributed  to  the  presence  of  metallic  contaminants  at  the  surface  of 
this  diode.  If  this  were  true,  the  DLTS  measurements  on  the  B-5  diode 
should  have  indicated  a higher  defect  density  than  was  found  on  the  B-8 
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diode.  DLTS  scans  performed  on  N+/Pb  diode  from  both  the  B-5  and  B-8 
die  indicate  the  presence  of  two  hole  traps  in  both  the  B-5  and  B-8 
N+/ Pb  BP  diodes.  Figure  8.11  shows  the  DLTS  scan  of  the  B-5  N+/ Pb  BP 
diode.  The  first  hole  trap,  denoted  as  H2,  has  its  energy  level  at  Ey  + 
0.47  eV.  The  second  hole  trap,  denoted  as  H4,  has  its  energy  level  at 
Ev  + 0.8  eV.  The  trap  density  of  the  H2  is  1014  cm-3,  and  is  about  the 
same  for  both  the  B-5  and  B-8  diodes.  H4  is  the  dominant  trap  in  both 
the  B-5  and  B-8  diodes.  The  density  of  the  H4  trap  in  the  B-5  diode  is 
5 x 1014  cm-3.  In  the  B-8  diode,  the  H4  trap  was  found  to  be  located 
near  the  surface  region  of  the  diode  at  a depth  of  less  than  0.7  urn. 

In  summary,  the  I-V,  C-V  and  DLTS  measurements  on  the  diodes  from 
the  B-5  and  B-8  dies  from  JFET2-1358  run  show  good  correlation  between 
reverse  leakage  current  and  the  density  of  the  corresponding  metallic 
contami nants . 

8.3.2  Results  and  discussion  on  phase  1 group  C (P1PC)  diodes 

I-V,  C-V  and  DLTS  measurements  were  performed  on  some  of  the  group 
C diodes,  and  the  results  are  shown  in  Figures  8.12  through  8.19.  The 
reverse  I-V  characteristics  curves  shown  in  Figure  8.13  indicate  that 
the  Nb/Pc  BP  diodes  have  consistently  higher  reverse  leakage  than  the 
Pjj/N-  BP  diodes  on  the  same  die.  The  high  reverse  leakage  current  is 
attributed  to  an  induced  n-channel  being  formed  in  the  lightly  doped  P 
collector. 

C-V  and  DLTS  measurements  were  performed  on  the  C-10  Nb/ Pc  BP  and 
Pb/N-  BP  diodes,  and  the  results  are  shown  in  Figures  8.14  through 
8.19.  Both  the  majority  carrier  (electron)  traps  in  the  N-  and  the 
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majority  carrier  (hole)  traps  in  the  Pc  were  analyzed  and  identified. 
In  the  C- 10  Pb/N-  BP  diode,  two  electron  traps  were  identified  (See 
Figs.  8.15-8.17).  The  first  electron  trap  El  has  an  energy  level  of 
Ec-0.19  eV  and  a density  of  1.7  x 1012  cm'3.  The  second  electron  trap 
E2  has  an  energy  level  of  Ec  - 0.56  eV  and  a density  of  9.3  x 1013 
cm  . Since  the  N-  dopant  density  is  one-tenth  of  the  Pb  dopant 
density,  the  observed  electron  traps  El  and  E2  should  be  in  the  N-  side 
of  the  junction. 

Two  shallow  hole  traps  were  observed  in  the  C-10  Nb/Pc  BP  diodes  as 
is  shown  in  Figure  8.18.  The  first  hole  trap  HI  has  an  energy  level  of 
Ev  + 0.22  eV  and  a trap  density  of  7.3  x 1013  cm'3.  The  energy  level  of 
the  second  hole  trap  H3  is  E y + 0.66  eV  (see  Figures  8.18  and  8.19). 
Table  8.4  shows  the  electron  and  hole  traps  observed  in  JFET  2-1358  run 
P1PB  and  P1PC  diodes.  A summary  of  the  deep-level  defects  observed  in 
P1PB  and  P1PC  diodes  is  given  in  Table  8.5. 

8.3.3  Results  and  discussion  on  phase  1 group  D (P1P2)  MOS 
capacitors 

The  results  of  the  C-V,  TSCAP  and  DLTS  measurements  on  group  D MOS 
capacitors  from  JFET2-1358  (die  #13,  #15,  and  #18)  run  are  presented  in 
this  section. 

C-V  measurements  vs.  temperature  (77  to  350  K)  and  frequency  (1  MHz 
to  50  MHz),  Thermal ly-Stimul  ated-Capacitance  (TSCAP)  and  Deep-Level - 
Transient-Spectroscopy  measurements  were  performed  on  the  phase  1 group 
D MOS  capacitors  with  N-substrate. 

The  results  are  quite  interesting,  and  conclusions  obtained  from 
this  study  are  consistent  with  the  results  obtained  from  the  DLTS  data 
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on  group  B and  group  C P-N  diodes  concerning  the  metallic  contaminants 
discussed  above.  The  results  of  our  analysis  on  the  C-V,  TSCAP  and  DLTS 
data  for  group  D MOS  capacitors  with  N-  substrate  are  summarized  as 
fol 1 ows : 

(1)  The  frequency  dependence  of  C-V  curves:  Fig.  8.20  shows  the 

C-V  plots  for  device  #18  for  four  different  frequencies 
(i.e.,  1,  5,  10,  20  MHz),  respectively.  It  is  noted  that  in 
the  accumulation  region  (i.e.,  V > 0)  the  insulator  (i.e., 
oxide)  capacitance  depends  strongly  on  frequency  (decrease 
with  increasing  frequency).  To  explain  this  result,  Fig. 

8.21  shows  the  equivalent  circuit  of  an  MOS  capacitor 
including  interface  trap  states  and  depletion  capacitance;  Cs 
and  Rs  are  associated  with  the  interface  trap  density.  The 
product  RSCS  is  defined  as  the  interface-trap  lifetime,  which 
determines  the  frequency  behavior  of  the  interface  traps. 
and  R.j  are  the  oxide  capacitance  and  oxide  shunt  resistance 
(which  is  infinity  for  ideal  case),  respectively.  The 
reduction  in  equivalent  capacitance  with  frequency  in  the 
accumulation  region  is  due  to  the  effect  of  R^ , R$  or  both. 

The  most  interesting  result  shown  in  Fig.  8.20  is  the 
hysteresis  loops  observed  in  the  inversion  region  of  the  C-V 
plots.  At  a reverse  bias  voltage  of  -20  V,  the  donor  states 
in  the  space  charge  region  become  positively  charged  due  to 
the  emission  of  electrons  from  the  deep  donor  traps  (as  will 
be  explained  later  in  the  DLTS  and  TSCAP  data),  and  increases 
drastically  to  a maximum  capacitance  determined  by  the  trap 
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density.  This  effect  is  more  pronounced  at  1 MHz  than  at 
other  higher  frequencies,  and  diminishes  at  20  MHz  in  which 
the  trap  states  can  no  longer  follow  the  ac  signal  (i.e,  RSCS 
time  constant  is  too  large).  The  results  also  showed  that 
the  capacitance  of  these  MOS  capacitors  becomes  constant  and 
independent  of  bias  voltage,  implying  that  the  interface 
traps  and  depletion  capacitance  are  negligible  at  very  high 
frequencies.  Fig.  8.22  shows  the  frequency  dependence  of  C-V 
curves  for  an  MOS  capacitor  on  p-col lector  substrate  at  300  K 
for  four  different  frequencies.  It  is  noted  that  in  contrast 
to  the  MOS  capacitor  built  on  the  n-  subtrate  shown  in  Fig. 
8.20  negligible  hysteresis  loops  were  observed  in  the  MOS 
capacitor  built  on  the  p-  collector  substrate. 

(2)  The  temperature  dependence  of  the  C-V  curves:  The  hysteresis 

loops  observed  in  the  inversion  region  of  the  C-V  curves  at 
300  K shown  in  Fig.  8.20  disappear  for  temperatures  below 
250K  (see  Fig.  8.23);  this  is  due  to  the  freeze-out  of  the 
donor  traps  (which  become  neutral  states).  Fig.  8.24  shows 
the  C-V  curves  vs.  frequency  (2  to  20  MHz)  for  capacitor  #18 
at  77  K.  It  is  clearly  shown  that  the  C-V  curves  are  quite 
normal  with  no  hysteresis  loops  observed  in  all  the  frequency 
ranges  studied.  Thus,  from  the  temperature  dependence  of  the 
C-V  characteristics  curves  it  is  obvious  that  the  hysteresis 
effect  observed  in  these  MOS  capacitors  is  due  to  the  donor 
traps  located  at  the  surface  of  the  N-substrate. 


123 


(3)  Thermal  ly-Sti mul ated-Capaci tance  (TSCAP)  measurements:  In 

order  to  explain  the  physical  reasons  for  causing  the 
observed  hysteresis  loops  shown  in  Fig.  8.20  and  Fig.  8.24, 
the  TSCAP  experiment  was  made  on  device  #18.  Fig.  8.25  shows 
the  TSCAP  scans  for  #18  MOS  capacitor  from  150  K to  350  K at 
1 and  3 MHz.  The  solid  lines  are  for  the  1 MHz  results,  and 
the  broken  lines  are  for  the  3 MHz  results.  The  curves  shown 
can  be  explained  as  follows:  If  the  capacitor  is  reverse 

biased  at  -20  V initially,  and  cooled  down  from  room 
temperature  to  liquid  nitrogen,  and  is  then  warmed  up  to  room 
temperature  under  the  same  reverse  bias  condition  as  the 
initial  condition,  we  should  see  the  result  (upper  solid 
curve)  shown  in  this  figure.  If  the  capacitor  is  initially 

positively  biased  so  that  all  traps  were  filled  by  electrons, 
then  the  TSCAP  scan  should  reveal  trap  information  from  a 
single  TSCAP  scan.  This  is  shown  in  the  lower  solid  curve  of 
the  C-V  plot.  A big  step  change  in  the  capacitance  was 
observed  in  device  #18  for  T > 250  K.  This  capacitance  step 
is  a result  of  the  electron  emisson  from  the  donor  trap 
states  at  the  S i - S i 0 2 interface  or  near  the  silicon 
surface.  The  broken  lines  are  the  TSCAP  scans  taken  on  the 
same  MOS  device  at  3 MHz  frequency.  Note  that  the  change  in 
capacitance  at  the  onset  of  the  thermal  emission  and  at  the 
end  of  electron  emission  is  directly  related  to  the  density 
of  the  electron  trap  via  the  relation: 


Nt  - (21C/C0)N„ 


(8.1) 


The  density  of  electron  trap  determined  by  the  above  formula 
and  AC  and  CQ  shown  in  Fig.  8.25  is  found  equal  to  3 x 10^ 
cm-3.  Fig.  8.26  shows  a similar  TSCAP  scan  of  the  same  MOS 
capacitor  performed  at  2 MHz  and  at  a reverse  bias  voltage  of 
-5  V.  It  is  noted  that  two  capacitance  steps  were  observed 
in  this  plot  between  250  and  280  K,  indicating  that  there  are 
two  electron  traps  located  very  closely  in  energy  in  this 
temperature  range.  This  is  verified  by  the  DLTS  scan  to  be 
discussed  next. 

(4)  The  DLTS  scan  of  electron  traps  was  also  made  in  #18  MOS 
capacitor,  and  the  result  was  shown  in  Fig.  8.27.  Two 
electron  traps  with  energies  of  0.4  and  0.56  eV  were  observed 
in  this  capacitor,  and  both  are  found  to  be  donor  type 
traps.  These  donor  type  traps  are  formed  due  to  the  metallic 
contamination  in  the  tested  devices. 

8.3.4  Results  and  discussion  on  phase  2 part  1 (P2P1)  diodes 
The  measured  results  of  D3  and  D4  diodes  in  phase  2 part  1 (P2P1) 
units  selected  from  run  JFET2-1358  are  shown  in  Table  8.6  and  8.7, 
Figures  8.28  through  8.31.  Table  8.6  and  8.7  summarize  the  experimental 
results  for  the  H3  (Ev  + 0.66  eV)  and  E2  (Ec  - 0.56  eV)  traps  observed 
in  the  tested  D3  and  D4  diodes.  The  trap  densities  for  the  H3  and  E2 
traps  are  in  the  range  of  4.1  x 1013  cm"3  and  1.9  x 1013  to  3.9  x 1013 
cm-3,  respectively.  Note  that  the  measured  trap  density  for  the  H3  and 
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E2  levels  was  lower  than  those  observed  in  the  units  reported  in  our 
phase  1 study.  The  reason  for  this  difference  has  been  attributed  to 
the  reserve  implant  procedure  used  in  these  two  units. 

The  reverse  I-V  curves  for  D3  diodes  of  unit  #14,  #19  and  #15 
selected  from  run  JFET2-1358  P2P1  wafer  lots  are  shown  in  Fig.  8.28,  in 
which  D3  diode  from  unit  #14  was  found  to  have  the  highest  leakage 
current  among  the  diodes  tested.  It  is  noted  that  #14,  03  diode  was 
also  found  to  have  the  highest  H3  trap  density.  Fig.  8.29  shows  the 
results  of  the  C-V  measurements  for  unit  #14  for  D3  and  D4  diodes, 
respectively,  for  frequencies  of  1 MHz  and  20  MHz.  The  results  of  C-V 
measurements  on  diodes  D3  and  D4  showed  that  the  capacitance  values  for 
03  depend  strongly  on  frequency  but  little  dependent  of  frequency  for  D4 
diode.  It  is  probably  because  the  D4  diode  has  low  E2  trap  density. 

The  DLTS  scan  of  main  hole  trap  H3  and  main  electron  trap  E2  in  the 
D3  and  D4  diodes  from  JFET2-1358  wafer  lots  are  also  shown  in  Figs.  8.30 
through  8.31,  respectively. 

In  order  to  find  the  correlation  between  the  metallic  contaminants 
induced  hole  trap  density  and  the  1/f  noise  in  seven  NPN  and  PNP 
transistors  as  well  as  JFETS  from  JFET2-1358  run,  the  1/f  noise  were 
measured  in  Harris  on  devices  in  which  DLTS  and  C-V  measurements  have 
also  been  made.  The  results  of  the  1/f  noise  measurements  at  five 
different  frequencies  (i.e,  10,  100,  IK,  10K,  and  100K)  on  seven  npn  and 
pnp  transistors  as  well  as  JFETs  are  listed  in  table  8.8,  8.9,  and 
8.10.  Note  that  the  highest  leaky  device  is  listed  at  the  top  of  table 
8.8,  8.9,  and  8.10,  while  the  lowest  leaky  device  is  listed  at  the 
bottom  of  the  table.  From  this  table  it  is  clearly  shown  that  there  is 
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good  correlation  between  the  hole  trap  density  (Ev  + 0.66  eV  hole  trap) 
and  the  leakage  current.  However,  it  should  be  noted  that  the  noise 
figures  for  both  the  NPN  transistors  and  the  JFETs  are  not  closely 
related  to  the  defect  density  deduced  from  the  DLTS  data.  The  only  good 
correlation  obtained  among  the  leakage  current,  trap  density,  and  the 
1/f  noise  figures  are  those  deduced  from  the  PNP  transistors,  as  is 
shown  in  Table  8.9.  In  order  to  find  a better  correlation  between  the 
noise  figures  and  trap  density  due  to  metallic  contamination,  it  is 
necessary  to  perform  the  noise  measurements  at  frequencies  below  10  Hz. 

8.4  Summary  and  conclusion 

From  the  analysis  of  the  I-V,  C-V,  TSCAP,  and  DLTS  data  on  the  P/N 
and  N/P  diodes  and  MOS  capactors  from  the  JFET2-1358  run,  we  reach  the 
following  conclusions: 

(1)  The  metallic  contaminants  presented  in  these  tested  devices 
introduce  several  deep  level  hole  and  electron  traps.  The 
main  hole  trap  is  due  to  the  Ev  + 0.66  eV  level  and  the  main 
electron  trap  is  due  to  the  Ec  - 0.56  eV  level  with  the 
average  trap  density  in  the  low  to  mid  1014  cm-3  range.  The 
surface  concentration  of  these  traps  is  expected  to  be  much 
higher  than  the  value  measured.  If  the  hole  trap  is  a donor 
type  defect,  it  will  have  a detrimental  effect  on  the  low- 
doped  p-type  surface  since  an  n-inversion  layer  may  be  formed 
by  a large  density  of  such  hole  trap  contaminated  at  the 
surface,  resulting  in  a large  surface  leakage  current  in  this 
n-channel  inversion  layer.  This  is  indeed  the  case  for  the 
lateral  N+/P  diodes  such  as  R-5  and  C-10  diodes  studied  here. 
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(2)  The  reverse  I-V  characteristics  were  found  closely  related  to 
the  DLTS  results  in  that  the  vertical  junction  diodes  showed 
very  little  leakage  current  and  smaller  density  of  metalllic 
contaminant  while  the  lateral  N+/P  diodes  showed  a very  large 
surface  leakage  current  and  a large  density  of  surface  defect 
density.  As  for  the  lateral  P+/N  diodes,  the  surface  leakage 
current  was  moderate  when  compared  to  the  lateral  N+/P 
diodes.  This  is  due  to  the  fact  that  N-surface  will  not  be 
inverted  by  the  donor-type  metallic  contaminants  while  an 
n-inversion  channel  may  be  formed  in  the  lightly  doped 
p-surface  layer  by  the  metallic  contaminant. 

(3)  The  results  of  our  C-V,  TSCAP,  and  DLTS  measurements  on  #18 
MOS  capacitor  are  quite  interesting  and  clearly  show  that  a 
majority  donor  type  electron  trap  located  at  the  Si -Si 02 
interface  plays  a key  role  in  influencing  the  temperature  and 
frequency  dependence  of  the  high  frequncy  C-V  behavior  in 
this  MOS  capacitor.  The  density  of  this  donor  trap  is  very 
large  (of  the  order  of  1 x 1015  cm-3),  and  is  believed  to 
have  a substantial  effect  on  the  1/f  noise  observed  in  the 
JFET  devices  fabricated  from  JFET2-1358  run.  The  donor  trap 
observed  in  this  MOS  capacitor  is  believed  to  be  due  to  the 
same  electron  trap  we  observed  in  the  lateral  P/N  diode 
(C-10)  discussed  earlier.  Thus,  this  donor  type  electron 
trap  may  be  related  to  the  metallic  contaminants  shown  in 
Table  8.5. 
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(4)  Since  all  the  electron  and  hole  traps  observed  in  the  diodes 
and  MOS  capacitors  from  implanted  JFET  wafer  lots  are  donor- 
like traps,  it  is  expected  that  p-silicon  surface  is  more 
likely  to  be  inverted  than  the  n-silicon  surface.  Therefore, 
an  n-type  inversion  channel  can  be  easily  formed  in  a lateral 
N+/P  diode,  resulting  in  a larger  leakage  current  and  1/f 
noise  in  such  a device. 

(5)  The  hysteresis  loops  observed  in  the  C-V  scans  for  some  of 
the  MOS  capacitors  from  the  implanted  JFET  wafer  lot  can  also 
be  attributed  to  the  trap  states  introduced  by  the  metallic 
contaminants  at  the  S i - S i 0 2 interface,  in  which  the  emission 
of  electrons  from  the  electron  traps  at  a characteristic 
temperature  results  in  the  hysteresis  behavior  observed  in 
the  C-V  plot. 

(6)  From  the  results  presented  in  section  8.3.4,  it  is  noted  that 
there  is  a good  correlation  between  the  trap  density  and  the 
reverse  leakage  current  observed  in  the  P2P1  devices, 
implying  that  the  reverse  leakage  current  observed  in  diodes 
assembled  from  JFET2-1358  run  was  due  to  the  bulk  generation 
recombination  current  stemmed  from  thermal  generation  and 
recombination  via  the  electron  and  hole  trap  centers  such  as 
E2  and  H3  traps  reported  above.  As  for  the  1/f  noise  source, 
the  results  for  the  P2P1  devices  shown  in  Table  8.8-8.10  are 
more  random  and  scattered  in  nature,  and  showed  only  a 
certain  degree  of  correlation  to  the  density  of  metallic 
contaminants  observed  in  these  devices. 
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In  short,  this  chapter  presents  the  results  of  our  studies  of  the 
effects  of  metallic  contamination  on  the  device  performance.  It  helps 
us  to  obtain  a better  understanding  concerning  the  correlation  among  the 
reverse  leakage  current,  1/f  noise  source,  and  the  density  of  metallic 
contaminants  determined  by  the  C-V  and  DLTS  measurements. 
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Table  8.1 


List  of  the  tested 

devices  used 

in  the  metallic  contamination  study. 

Group  B 

D1  (N+/Pb)  diode 
D2  (P+/Nb)  diode 

B-5,  B-8 

Part  1 

(PIPb) 

Group  C 

D3  ( N^/ Pq ) diode 

C-9,  C-10 

Phase  1 

(P1PC) 

D4  (Pb/N_)  diode 

Group  D 

Part  2 

( P1P2) 

MOS  Capacitor 

#13,  #15,  #18 

D3  (Nb/Pc)  diode 

#6,  #11,  #14, 
#15,  #18,  #19 

Phase  2 Part  1 

(P2P1) 

#20,  #28,  #31, 

D4  (Pb/N_)  diode 

#32,  #34,  #38, 
#39. 
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Table  8.2 


Results  of  the  I-V  and  C-V  measurements  on  Group  B diodes. 


DIE  PARAMETERS 


N+/Pb  BP  P+/Nb  BP  N+/Pb  SW 


B-5  Diode  Ideality  Factor 
Base  Concentration 


1.05  1.06 

2 x 1016  9.7  x 1016 


1.06 

2.6  x 1017 


B-8 


Diode  Ideality  Factor 
Base  Concentration 


1.10  1.03  1.04 

1.8  x 1016  1.2  x 1017  2.5  x 1017 
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Table  8.3 


Results  of  the  I-V  and  C-V  measurements  on  Group  C diodes. 


DIE  PARAMETERS 


Pb/N-  BP  Nb/Pc  BP  Nb/Pc  SW 


C-9  Diode  Ideality  Factor  1.12  1.10  1.06 

Collector  Concentration  2.1  x 10^  1.1  x 10^  1.3  x 10^ 


C-10  Diode  Ideality  Factor  1.10  1.10  1.04 

Collector  Concentration  2.0  x 10*5  1.2  x 10^  1.0  x 10^ 


133 


Table  8.4 


Electron  and  hole  traps  observed  in  P1PB 
and  P1PC  diodes  from  JFET2-1358  run. 


DIODES 

ELECTRON 

TRAPS 

HOLE 

TRAPS 

Ec  - et 

NTo 

EV  ' 

- et 

"T, 

(eV) 

(cm  3) 

(eV) 

(cm  ) 

B-5  N+/Pb 

BP 

H2  = 

0.47 

1.0 

x 1014 

H3  = 

0.66 

1.0 

x 1014 

H4  = 

0.83 

5.0 

x 1014 

B-8  N+/Pb 

BP 

H2  = 

0.47 

1.0 

x 1014 

H3  = 

0.66 

1.0 

x 1014 

H4  = 

0.83 

1.0 

x 1014 

C-10  Pb/N. 

. BP 

El  = 0.19 

1.7  x 1012 

E2  = 0.56 

1.0  x 1014 

CIO  Nb/Pc 

BP 

HI  = 

0.22 

7.3 

x 1013 

H3  = 

0.66 

5.3 

x 1014 

H4  = 

0.83 

8.9 

x 1013 
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Table  8.5 


Summary  of  deep-level  defects  observed  in  the  P1PB  and  P1PC  diodes. 


JFET  DEVICES 

DEFECTS 

ET(eV) 

DLTS  PEAK  (K) 

en  or 

ep  (s_1) 

°n  or 
op  (cm2) 

Fe 

El  = Ec- 

0.19 

138 

555.6 

7 x 10-17 

(Au,  Co,  Ti) 

E2  = Ec- 

0.56 

295 

555.6 

1.2  x 10"16 

Fe 

HI  = Ev+ 

0.22 

152 

86 

1.4  x 10"16 

Fe 

H2  = Ev  + 

0.47 

217 

277 

9 x 10“13 

H3  = Ev+ 

0.66 

334 

86 

1.2  x 10'14 

H4  = Ev  + 

0.8  ~ 

0.86  369 

86 

4.7  x 10-13 
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Table  8.6 

Defect  density  and  background  dopent  density  for  the 
N^/Pc  and  P^,/ Nc  diodes  from  JFET2-1358  processes. 


*Dev i ce 

Parameters 

Background  Density 
(cm"3) 

Trap  Density 
(cm-3) 

#14 

(red) 

(high) 

Vpc 

VNc 

1.03  x 1016 
2.7  x 1015 

$H3  = 1.07  x 1014 
&E2  = 2.2  x 1013 

#15 

(red) 

(low) 

Vpc 

VNC 

1.23  x 1016 
3.0  x 1015 

H3  = 4.8  x 1013 
E2  = 1.9  x 1013 

#19 
(bl ue) 
(medium) 

Vpc 

VNC 

1.03  x 1016 
3.1  x 1015 

H3  = 6 x 1013 
E2  = 2.2  x 1013 

#32 
(bl ue) 
(high) 

Vpc 

VNC 

9.8  x 1016 
2.86  x 1015 

H3  = 9.13  x 1013 
E2  = 3.3  x 1013 

#38 

(green) 

(medium) 

Vpc 

1.07  x 1016 
3.0  x 1015 

H3  = 8.1  x 1013 
E2  = 3.9  x 1013 

#39 

(green) 

(high) 

Vpc 

VNC 

1.09  x 1016 
2.7  x 1015 

H3  = 1.03  x 1014 
E2  = 3.4  x 1013 

* Area:  1.652  x 10-3  cm3 
$ H3:  Ev  + 0.66  eV 

& E2:  Ec  - 0.56  eV 
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4.48  1.96  0.88  0.52  1.19 


Hole  trap  density,  leakage  of  current  and  the  1/f  noise  data  for  P2P1  PNP  transistors  from  JFET2-1358  run. 
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Table  8.10 


Hole  trap  density,  leakage  current  and  1/f  noise  data 
for  P2P1  JFETs  from  JFET2-1358  run. 


SN 

Leak  rank 
(Hi  to  Lo) 

Hole  density 
(cm-3) 

10 

JFET 

100 

en  (nV/(Hz)1/2) 
IK  10K 

100K  (Hz] 

#18 

47.9 

12  x 1013 

140 

85 

50 

20 

10 

#28 

48.1 

11  x 1013 

165 

100 

60 

22 

11 

#20 

102 

8.4  x 1013 

205 

115 

64 

23 

11 

#11 

104 

7.2  x 1013 

130 

82 

49 

21 

10 

#34 
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5.4  x 1013 

190 

100 

59 

21 

10 

#31 
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4.7  x 1013 
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64 

21 

10 

#6 

145 

4.1  x 1013 

140 

99 

57 

22 

11 

140 


Figure  8.1  Photolithographic  masks  used  in  the  fabrication  of 

test  devices  for  the  I-V,  C-V,  and  DLTS  measurements. 


Structure  Dimensions:  N-base  16  mils  x 16  mils 


Figure  8.2  Test  structure  for  a Nb/Pc  bottom  plate  diode(BP). 
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Structure  Dimensions:  N-base  stripe  16  mils  x 0.6  mils 


Figure  8.3  Test  structure  for  the  N^/Pc  sidewall  diode(SW). 
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D3  (Nb/Pc) 


Figure  8.4  Cross  sectional  views  of  the  Dl,  D2,  D3,  and  D4 
diodes. 
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POLY  POLY 


( Not  To  The  Scale  ) 

Figure  8.5  Top  and  cross  sectional  views  of  a MOS  capacitor  from 
#JFET  2-1358  run. 
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Figure  8.6  Forward  current-voltage  curves  for  three  B-8  diodes 
selected  from  JFET  2-1358  (P1PB)  run. 
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Figure  8.7  Forward  current-voltage  curves  for  three  B-5  diodes 
selected  from  #JFET  2-1358  (P1PB)  run. 
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Figure  8.8  Reverse  current-voltage  curves  for  B-5  and  B-8 
diodes  from  #JFET  2-1358  (P1PB)  run. 
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1/C2  vs.  VR  for  B-5  and  B-8  diodes  from  #JFET 
2-1358  (P1PB)  run. 


Figure  8.9 
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Vr(V) 


Figure  8.10  1/C2  vs.  VD  for  group-B  and  -C  diodes  selected  from 

#JFET  2-1358  (P1PB  and  P1PC)  run. 


DLTS  SIGNAL 
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Figure  8.11  DLTS  scan  of  hole  traps  in  the  B-5  N+/P.  BP  diode 
selected  from  #JFET  2-1358  (P1PB)  run. 
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vF(v) 


Figure  8.12  Forward  I-V  curves  for  C- 10  diodes  selected  from 
#JFET  2-1358  (P1PB)  run. 
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Figure  8.13  Reverse  I-V  curves  for  group-C  diodes  selected  from 
#JFET  2-1358  (P1PC)  run. 
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Fiaure  8 14  Capacitance-voltage  curves  for  C-10  diodes  measured 
at  20  MHz  from  #JFET  2-1358  (P1PC)  run. 
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TEMPERATURE  (K) 


Figure  8.15  DLTS  scan  of  electron  traps  in  the  C-10  P./N-  BP 
diode  selected  from  #JFET  2-1358  (P1PC)  run. 
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Figure  8.16  Thermal  electron  emission  rates  divided  by  T vs. 

inverse  of  temperature  for  El  electron  trap  in  the 
C-10  P./N-  BP  diode  selected  from  #JFET  2-1358 
( P1PC)  run . 


Figure  8.17  Thermal  electron  emission  rates  divided  by  T'1  vs.  inverse  of  temperature  for  E2 
electron  trap  in  the  C-10  P5/N-  diode  selected  from  #JFET  2-1358  (P1PC)  run. 


156 


en/T 


2 


) 


DLTS 


157 


100  150  200  250  300  350 

TEMPERATURE  (K) 


Figure  8.18  DLTS  scan  of  hole  traps  in  the  C-10  N./P  BP  diode  selected 
from  #JFET  2-1358  (P1PC)  run. 
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2 

Figure  8.19  Thermal  hole  emission  rates  divided  by  T vs.  inverse 
of  temperature  for  H4  hole  trap  in  the  C-10  N./P  BP 
diode  selected  from  #JFET  2-1358  (P1PC)  run. 
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Figure  8.20  C-V  curves  vs.  frequency  at  T = 300  K for  #18  MOS  ( N- 
substrate  ) capacitor  selected  from  run  #JFET  2-1358  ( 
P1P2  ).  Hysteresis  loop  is  clearly  shown  for  f < 10  MHz 
and  V < -5 V. 
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T - cs  *s 


Cp  = CD  + Cs/(  1 +,j2  ^ 1 

G = C w2t/(  1 + co2  x2  ) 

P * 


(a) 


(b) 


Figure  8.21  Equivalent  circuit  including  interface-trap  effect, 
where  C and  R are  associated  with  interface  trap 
density! 
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Figure  8.22  C-V  curves  vs.  frequency  at  T = 300  k for  #15  MOS 

capacitor  on  P -substrate  selected  from  #JFET  2-1358 
(P1P2)  run.  Hysteresis  effects  is  negligible  in 
this  capacitor. 
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Figure  8.23  C-V  plots  at  1 MHz  as  a function  of  temperature  for 
#18  MOS  capacitor  with  N-substrate  selected  from 
#JFET  2-1358  (P1P2)  run.  Note  that  hysteresis 
effect  disappears  for  $ 250  k. 
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Figure  8.24  C-V  plots  at  77  K for  #18  MOS  capacitor  (N-substrate) 
selected  from  run  #JFET  2- 1 358 ( PI P2 ) for  f = 2,  5,  10, 
and  20  MHz.  Note  that  no  hysteresis  loop  is  shown  in 
the  C-V  curves. 
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Figure  8.25  The  thermally-Stimul ated-Capaci tance  (TSCAP)  scans  for 
#18  MOS  capacitor  (N-substrate)  at  1 and  3 MHz  from  run 
#JFET  2-1 358 ( PI P2 ) . The  capacitance  step  observed  for 
T > 265  K is  due  to  the  electron  emission  from  the  trap 
in  the  device. 
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Figure  8.26  TSCAP  scan  taken  at  2 MHz  for  #13  MOS  capacitor  (N- 

substrate)  showing  electron  traps  for  250  < T < 300  k. 
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Figure  8. 27 


DLTS  scan  of  electron  traps  for  #18  MOS  capacitor  with 
N-substrate  selected  from  run  #JFET  2-1 358 (PI P2) . Two 
donor-type  electron  traps  were  detected. 


167 


5. 200/di v (V) 
VR 


Figure  8.28  Reverse  I-V  curves  for  D3  diodes  of  unit  #14,  #19  and  #15 
selected  from  run  #JFET  2-1358  (P2P1)  wafer  lots. 
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Figure  8.29  Capacitance  vs.  VR  for  diodes  D3  and  D4  of  unit  #14 
from  #JFET  2-1358 ( P2P1 ) run  wafer  lot. 
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Figure  8.30  DLTS  scan  of  hole  trap  H3  in  the  unit  #39  D3  diode 
selected  from  #JFET  2-1358  (P2P1)  run  wafer  lot. 
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Figure  8.31  DLTS  scan  of  the  electron  trap  E2  in  the  unit  of 
#32  D4  diode  selected  from  #JFET  2-1358  (P2P1) 
run  wafer  lot. 


CHAPTER  IX 

SUMMARY  AND  CONCLUSION 


A theoretical  study  of  the  effects  of  increasing  the  density  of 
shallow  impurity  dopants  in  a degenerate  semiconductor  device  and  an 
experimental  study  of  the  metallic  contamination  on  device 
characteristics  has  been  presented  in  this  dissertation.  The  highlights 
and  conclusions  from  this  study  are  summarized  as  follows: 

(1)  Expressions  for  the  Debye  length,  built-in  potential  and 
depletion  layer  width  of  a heavily-doped  pn  junction  diode  as 
well  as  a heterojunction  diode  are  reformulated  by  including 
the  heavy  doping  effects  such  as  carrier  degeneracy,  bandgap 
narrowing,  and  dopant  density-dependent  dielectric 
constant.  The  results  from  numerical  calculations  show  that 
a)  bandgap  narrowing  and  carrier  degeneracy  have  important 
effects  on  the  junction  built-in  potential;  b)  the  dopant 
impurity  density-dependent  constant  is  the  only  heavy  doping 
effect  which  will  have  a strong  influence  on  the  depletion 
layer  width  and  depletion  capacitance;  c)  the  effect  of 
carrier  degeneracy  will  be  the  dominant  heavy  doping  effect 
in  the  heavily-doped  GaAs-AlGaAs  devices  due  to  the  very 
small  electron  effective  masses  in  GaAs  and  AlGaAs. 

(2)  The  generalized  Einstein  relation  for  the  majority  carriers 
and  minority  carriers  in  a heavily-doped  semiconductor  are 
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still  valid  when  assuming  the  parabolic  density  of  states  and 
rigid  band  approximations . Additional  terms  should  be  added 
in  when  considering  the  impurity  band  broadening  and  band 
tail  states  instead  of  employing  the  parabolic  density  of 
states  and  rigid  band  approximation. 

(3)  A closed  form  expression  for  the  I-V  characteristics  of  a 

heavily-doped  Si  n+p  graded  junction  is  presented  by 

modifying  the  theory  of  anisotropic  heterojunctions , in 
which,  a transition  region  is  modeled  to  extend  beyond  the 
space  charge  region  with  a costant  bandgap  gradient. 
Numerical  analysis  show  that  the  effects  of  heavy  doping, 
especially  the  bandgap  narrowing  and  the  carrier  degeneracy, 
have  a significant  influence  on  the  electron  current  density 
at  the  quasi -neutral  region  of  a heavily-doped  pn  graded 
junction.  This  is  because  the  bandgap  gradient  affects  the 
transport  of  the  injected  minority  carriers  by  influencing 
the  carrier  distribution  in  the  quasi -neutral  regions. 

(4)  The  theoretical  analysis  of  the  effects  of  heavy  doping  on 
the  electrical  characteristics  for  GaAs  JFETs  and  MESFETs 
show  that  the  carrier  degeneracy  has  a significant  effect  on 
the  threshold  voltage  and  the  I-V  characteristics  of  heavily- 
doped  E-JFETs  and  E-MESFETs,  while  the  bandgap  narrowing  has 
negligible  influence  on  the  electrical  characteristics 
because  of  the  very  small  electron  effective  mass. 

(5)  From  the  analysis  of  I-V,  C-V,  TSCAP,  and  DLTS  data  on  the 
metallic  contaminated  silicon  diodes  and  MOS  capacitors,  we 
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conclude  that  a)  the  metallic  contaminants  presented  in  the 
tested  devices  will  introduce  several  deep  level  hole  and 

electron  traps.  The  main  hole  trap  is  due  to  the  Ev  + 0.66 

eV  level  and  the  main  electron  trap  is  due  to  the  Ec  - 0.56 

eV  level  with  the  average  strap  density  in  the  low  to  mid 
1014  cm-3  range;  b)  the  results  of  I-V  and  DLTS  measurements 
show  that  the  reverse  leakage  current  is  closely  related  to 
the  density  of  metallic  contaminant  induced  traps  in  the 
corresponding  tested  devices;  c)  experimental  evidence  from 
the  C-V,  TSCAP,  and  DLTS  measurements  indicate  that  the 
metallic  contamination  will  influence  the  temperature  and 
frequency  dependence  of  the  high  frequency  C-V  behavior  in 
the  tested  MOS  capacitors.  An  abnormal  hysteresis  behavior 
was  also  observed;  d)  the  results  of  the  measured  1/f  figures 
showed  only  a certain  degree  of  correlation  to  the  density  of 
metallic  contaminant  induced  deep  level  trap. 

The  theoretical  and  experimental  studies  presented  in  this 
dissertation  will  certainly  help  us  to  obtain  a better  understanding  of 
the  effects  of  heavy  doping  and  metallic  contamination  on  the  electrical 
characteristics  of  silicon  devices,  and  the  theoretical  reults  derived 
from  this  research  may  be  applied  to  study  other  III -V  compound 


semiconductor  devices. 


APPENDIX 

DERIVATION  OF  THE  DEPLETION  LAYER  WIDTH  VS.  BUILT-IN 
POTENTIAL  FOR  A HEAVILY  DOPED  LINEARLY  GRADED  PN  JUNCTION 


To  derive  an  expression  for  the  depletion  layer  width,  x^,  vs  the 
modified  built-in  potential,  4^,  in  a heavily  doped  linearly  graded  pn 
junction  taking  into  account  of  the  effects  of  carrier  degeneracy, 
impurity  density-dependent  dielectric  constant,  and  bandgap  narrowing, 
we  consider  a heavily-doped  linearly  graded  junction  with  the  space 
charge  regions  located  at  _xpi_  an<^  xn|_»  anc*  a ne^  dopant  concentration 
gradient,  a.  Thus,  the  net  dopant  density  in  this  linearly  graded 
junction  may  be  represented  by  Np  - NA  = ax.  From  discussions  given  in 
section  2.1,  the  dopant  impurity  density-dependent  dielectric  constant, 
K$n(x)  and  Ksp(x),  for  a linearly  graded  pn  junction  can  be  expressed  by 


K (x)  = A e 
snv  ' n 


B N 
n xn 


0 < x < x 


H 

nL 


(A-l) 


00 


B N 
, P XP 


< x < 0 


( A— 2 ) 


where,  Nxn  = ND  - NA  = ax,  and  Nxp  = NA  - N0  = -ax.  The  Gauss  law, 
including  the  effect  of  impurity  density-dependent  dielectric  constant, 
is  given  by 


V (Ks(x)  E(x) ) = p/eQ 


(A-3) 
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Since  E(x^)  = ^(xpb)  = we  obtain 


qa  1 f H2  2 \ 
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(A-4) 
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Thus,  the  potential  can  be  expressed  by 
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*(x)  = ♦n  " 2^  ^ nL  A exp(B  'ax)  (x 
ox  n n ' 
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x2)  dx  ; 0 < x < x^  (A-6) 
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*(x)  - *p  - 2®.  j Pp  exp(-Bp  ax) 


(x2  ' xpL^  dx  ’ "xpL  x x x 0 


where,  we  have  substituted  Eqs.  (A-l)  and  (A-2)  into  Eqs.  (A-4)  and 
(A-5).  By  using  integration  by  parts,  Eqs.  (A-6)  and  (A-7)  are  reduced 
to : 
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where,  Cn  = -Bna,  and  Cp  = Bpa.  Thus,  the  built-in  potential,  4>bi , tor 
the  heavily-doped  linearly  graded  pn  junction  can  be  obtained  from  Eqs. 
(A-8)  and  (A-9)  as 
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where  determination  of  xdL  is  made  possible  by  iterating  Eqs.  (12)  and 
(A-10).  It  is  noted  that  values  of  K$n  and  Ksp,  which  appeared  in  Eq. 
(A-10)  in  terms  of  An>  Bn,  Ap,  and  Bp,  are  determined  from  fitting  the 
experimental  data  [2]. 
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